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1. Introduction 

1.1 Motivation 

One of the main motivations for studying integrable lattice models is their role as a lattice regular- 
ization of quantum field theories in continuous space-time. Integrable nonlinear sigma models are 
of particular interest, and recently there has been a growing interest in nonlinear sigma models with 
non-compact target spaces. This interest is motivated by possible applications to string theory on 
curved space-times in general, and to gauge theories via the AdS-CFT correspondence in particular. 

However, the quantization and the solution of such non-compact nonlinear sigma models still rep- 
resents a major challenge for the field of integrable models. Compared to the better understood 
nonlinear sigma models with compact target spaces one may expect important qualitative differ- 
ences, which make it problematic to apply the known techniques from the compact cases to the 
sigma models with non-compact target spaces. This point is exemplified by the relation between the 
Wess-Zumino-Novikov-Witten (WZNW) models associated to compact and non-compact sym- 
metric spaces respectively. The solution of the latter is possible ITU , but it is considerably more 
difficult than the solution of WZNW models associated to compact groups. 

In more general sigma models one can not hope to find the powerful Kac-Moody symmetries of 
the WZNW models but the integrable structure may still survive. In order to enter the next level of 
complexity one may therefore try to exploit the integrability of some of these models. Turning to 
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mental Research grant 05-01-00922, and a grant from the Russian Science Support Foundation 
2 Part of this work was supported by a DFG Heisenberg fellowship 
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a new class of models it is always advisable to look for the simplest member which still exhibits 
most of the new qualitative features. In the case of the conformal WZNW models it has turned out 
that Liouville theory already displays many of the relevant differences which distinguish the non- 
compact WZNW models from rational conformal field theories |T2|. Moving outside of the class 
of sigma models soluble thanks to Kac-Moody or similarly powerful chiral symmetries it seems 
natural to look for a useful counterpart of Liouville theory within this larger class of models. 

A natural candidate for such a model exists: the sinh-Gordon model. Indeed, there is some evidence 
BZZI ILul that the sinh-Gordon model can be seen as a "deformation" of Liouville theory which 
preserves its integrable structure when the conformal symmetry is lost. While there certainly exists 
a good basis for the study of the sinh-Gordon model in infinite volume — S-matrix and the form 
factors are known IIVGI IFMS , KMul iBLl ILel and the basic ingredients of the QISM approach were 
developed IS II — there does not seem to exist a sytematic approach to the quantization and solution 
of the sinh-Gordon model in finite spatial volume yet. Part of the problem is due to the usual 
divergencies and ordering problems of quantum (field) theory. But the other part of the problem 
seems to be closely related to the non-compactness of the target space in the sinh-Gordon model. 

Our main motivation behind the present project was therefore to find an integrable lattice regulariza- 
tion for the sinh-Gordon model. This not only tames the usual short distance singularities, it will also 
allow us to take care of the troubles from non-compactness of the target space in a mathematically 
well-defined framework. One particular feature that directly follows from the non-compactness of 
the target space will be the failure of the usual algebraic Bethe ansatz method [Fl | for the model at 
hand. This failure means that we will have to use the more general separation of variables method 
ltS2llS3l [SmT1 instead. 

1.2 Lattice sinh-Gordon and the modular XXZ magnet 

A quantum integrable system is a quantum system (TL , A, H), with Hilbert space H, algebra of ob- 
servables A, Hamiltonian H, in which there exists a set Q = {To, Ti, . . . } of self-adjoint operators 
such that 



Property (C) expresses completeness of the set Q of integrals of motion. It is equivalent to the 
statement that the spectrum of Q is non-degenerate, i.e., that simultaneous eigenstates of T&, k G 
Z-° are uniquely determined by the tuple of their eigenvalues. We will consider the so-called one- 
dimensional lattice models for which one has 



(A) 
(B) 
(C) 



[T,T']=0 VTJ'eQ, 
[T, H] = V T € Q, 

if [T, 0] = for all TeQ, then = 0(Q). 



H = K. m , 



A = B m , 



(1) 



with one copy of Hilbert space K and algebra of local observables B being associated to each of the 
N sites of a one-dimensional lattice. 
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The quantum inverse scattering method (QISM) HFSTl fFl I goes a long way towards the construction 
of large classes of quantum integrable models of this type. In this framework one usually charac- 
terizes JC as a representation of a Hopf algebra it of "symmetries", and B is generated from the 
operators which represent the elements of il on JC. It is clear that the representation theoretic prop- 
erties of JC will influence the physical properties of the resulting integrable model decisively. Good 
control over these properties will be crucial in the construction and solution of such models. 

In general it is a highly nontrivial problem to find the "right" representation JC which leads to a useful 
lattice regularization of a particular quantum field theory. We will here propose a particular choice 
for JC which will lead to a lattice model with particularly nice properties, and which will be shown 
to yield the sinh-Gordon Hamiltonian density in the continuum limit of the corresponding classical 
lattice model. The representations in question will be representations of the non-compact real form 
U q (sl(2,R)) of U q {sh) which have been studied in llPTTI lF3l |PT2l [BTl . 

The non-compactness of the target space will be reflected in the infinite-dimensionality of the repre- 
sentation JC. It is furthermore worth noting that the same representations were previously found to re- 
flect a key internal structure of Liouville theory IIPT1IIT2I . In view of the existing evidence MZZI ILul 
for the connection between Liouville theory and the sinh-Gordon model, it is quite natural that the 
same class of representations appears in our lattice version of the sinh-Gordon model as well. 

The corresponding representations possess a remarkable duality — they are simultaneously repre- 
sentations of U q {s\2) and ^(sfe), where q — e lnb and q = e Z7rb . One may therefore view them 
llF3l as representations of the modular double U q (sl2)®l{q(sl2) (see also BKLSl FBTI). The parameter 
b turns out to be proportional to the coupling constant (3 of the sinh-Gordon model. The self-duality 
of our representations will be directly related to the self-duality of the sinh-Gordon model under 
b — > b^ 1 which was previously observed in its scattering theory. The importance of this self-duality 
for our analysis can hardly be over-emphasized. 

It turns out that there is a close relative of our lattice sinh-Gordon model which is simpler in some 
respects. This integrable lattice model can be seen as a non-compact counterpart of the XXZ model 
with spins in infinite-dimensional representations of the modular double. We will refer to this model 
as the modular XXZ magnet. As some technical issues are simpler in the case of the modular XXZ 
magnet, we will first construct the latter model before we turn to the lattice sinh-Gordon model. 
In any case, it seems to us that the study of the modular XXZ magnet is of interest in its own 
right. We note in particular that despite the different underlying representation theory, our model has 
many structural similarities with the non-compact XXX type magnet based on infinite-dimensional 
highest weight representations of which was studied in [DKM KM|. The latter model plays an 
important role in high energy QCD JLTJ IFK1I . 

1.3 Plan of the paper 

To make our paper accessible to a reasonably wide audience, we presented the general description 
of our approach, the main definitions and results in the main body of the paper and collected more 
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technical developments in Appendices. The article is organized as follows. 

In Section 2 we define the modular XXZ magnet in terms of the representations V s describe its 
Hilbert space of states, construct the corresponding fundamental R-operator R(u), and discuss the 
construction of its Hamiltonian and the set of integrals of motion Q. 

The same is done for the lattice sinh-Gordon model in Section 3. We show that the Hamiltonian 
density of the sinh-Gordon model is recovered in the continuum limit of the corresponding classical 
model. We also show that the algebraic Bethe ansatz fails due to the non-compactness of the target 
space in the sinh-Gordon model. 

An important first step towards the solution of these models is taken in Section 4. We refine the 
spectral problem for the integrals of motion T& by constructing the Q-operator Q(u) which is related 
to the Tfc via the so-called Baxter equation. Analyzing the properties of Q(u), we derive a set 
of conditions for its eigenvalues qt(u) which can be seen as quantization conditions and which 
replace the usual Bethe ansatz equations in our models. The self-duality of our representations 
furthermore allows us to derive the so-called quantum Wronskian relation for the sinh-Gordon 
model with odd N, which encodes valuable additional information about the spectrum. 

In order to show that the conditions found in Section 4 are also sufficient to characterize the spectrum 
we apply the separation of variables approach to our models in Section 5. 

Section 6 contains concluding remarks on the conditions which characterize the spectrum of our 
models, the continuum limit, and the relation with the lattice and continuum versions of Liouville 
theory. We observe in particular that our results are consistent with the results and conjectures of 
BZallLul on the continuum sinh-Gordon model in a nontrivial way. 

Appendices contain necessary technical details. Appendix A collects the relevant information on the 
special functions that we use. Appendix B discusses the precise mathematical nature of the self- 
duality b — > b^ 1 of the representations that we use. Appendic C contains some important technical 
results on the structure of the monodromy matrix. Appendix D is devoted to the construction of the 
fundamental R-operator, the key object for the construction of local integrals of motion for the lattice 
models. Appendix E contains details on the derivation of the properties of the Q-operator Q(u). 

Acknowledgements. We thank S.Derkachov, L.Faddeev, S.Lukyanov, N. Reshetikhin, and especially 
F. Smirnov for stimulating discussions. A.B. is grateful to J.T. and R. Schrader for hospitality during his visits 
to the Institute for Theoretical Physics, FU-Berlin and to V. Schomerus for hospitality during a visit to DESY, 
Hamburg. 



5 



Contents 

1 Introduction [l] 

1.1 Motivation [l] 

1.2 Lattice sinh-Gordon and the modular XXZ magnet 

1.3 Plan of the paper 

2 Modular XXZ magnet Q 

2.1 Quantum group symmetry U q (sl(2, R)) 

2.2 Representations — algebra of observables 

2.3 Integrals of motion § 

2.4 Fundamental R-operator and Hamiltonian § 

2.5 Classical limit [Tl| 

2.6 Comparison with similar models [3 

3 Lattice sinh-Gordon model [3 

3.1 Definition of the model [3 

3.2 Fundamental R-operator and Hamiltonian [3 

3.3 Relation with the continuum theory 

3.4 Failure of the algebraic Bethe ansatz []3 

4 Q-operator and Baxter equation [l<] 

4.1 Explicit form of Q(w) [H 

4.2 Analytic properties of eigenvalues of Q (it) 

4.3 Self-duality and quantum Wronskian relation [2j] 

4.4 Parity and cyclic shift 

5 Separation of variables 

5.1 Operator zeros of B(m) 

5.2 Operators A(u) and D(u) 

5.3 Sklyanin measure 13 

5.4 Remaining cases [3 

6 Concluding remarks — outlook [5^ 

6. 1 On the Baxter equations 

6.2 Continuum limit 



6 

6.3 Connection with lattice Liouville model E3 

A Special functions E^l 

A. 1 Double Gamma function E3 

A.2 Function Wb(x) E3 

A.3 Function D a (x) El 

A.4 Integral identities for D a (a;) 

B Positivity versus self-duality of the representations V s S3 

C Structure of the monodromy matrix 

C. 1 Expansions in the spectral parameter 

C. 2 Quantum determinant 

D Construction of the fundamental R-operator R(u) 53 

D. 1 Clebsch-Gordan maps 

D.2 R-operator and braiding 

D.3 Proof of Theorem^ El 

D.4 Clebsch-Gordan and Racah-Wigner coefficients for V s E3 

D. 5 Proof of Proposition!] El 

E Construction of the Q-operator Q(u) E3 

E. 1 Preliminaries E3 

E.2 Construction of Q-operators El 

E.3 Proof of commutativity relations E2 



7 



2. Modular XXZ magnet 

In this section we will begin to develop the QISM for the modular XXZ magnet — an XXZ type 
non-compact spin chain, which has U q (sl(2, R)) as a quantum symmetry. 

2.1 Quantum group symmetry li q {si{2, R)) 

Let q = e* 7 , 7 = nb 2 , b € (0, 1). We will also use the notation Q = b + b~ l . 

The quantum group U q (s^) is a Hopf algebra with generators E, F, K, K^ 1 satisfying the relations 

KE — qEK, KF = q~ 1 FK, [E, F] — -^^{K 2 — K~ 2 ) (2) 

and equipped with the following co-product: 

ME) = E®K + K- 1 ®E, 

, A(K)=K®K. (3) 

A(F) =F®K + K~ 1 ®F, 

The relevant real form of U q (s\.2) is U q (sl(2, R)), which is defined by the following star-structure: 

K* = K, E* = E, F* = F. (4) 
The center of U q (sl(2, R)) is generated by the g-Casimir element: 

C= (2sm~f) 2 FE - qK 2 - q^K- 2 + 2, C* = C . (5) 

2.2 Representations V s — algebra of observables 

A one-parameter family of unitary representations V s of U q (5 1(2, R) ) can be constructed from a pair 
of self-adjoint operators p and x on L 2 (R) which satisfy [p, x] = (2m)~ 1 as follows: 

n s (E) ee E s =e+^ C ° Sh7rb(P ~ S) e^ b \ 



v,T7o. ^ = e ^ bp - (6) 

« S (F) ee F s =e -, bx cosh^(P + g ) e _, bX; 

sin 7 

For this representation 

C s ee tt s (C) = 4cosh 2 7r6s. (7) 

It is remarkable and important that the operators E s , F s and K s are positive self-adjoint. Indeed, the 
representations Vs are the only "reasonable" representations of U q (sl(2, R)) which have this prop- 
erty. This property will play a key role in much of the following developments. It will in particular 
ensure seld-adjointness of operators such as the Hamiltonian and the integrals of motion. It is also 
the mathematical basis for the self-duality of the representations V s , as shown in Appendix lBl(see. 
in particular, eq. dB.5b ). 
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The lattice model that we are about to define will have one of the representations V s attached to each 
site of the one-dimensional lattice. This means that we take 

H = {L 2 (R)) m (8) 

as the Hilbert space of our model, and let 

A=(ir s {U)) m i U=U q {sl{2,R)) (9) 

be a set of generators for our algebra of observables. Note that the operators in A are all unbounded, 
but there exists a basis for A whose elements are positive self-adjoint (see AppendixlBk The latter 
fact allows us to construct large classes of non-polynomial operator functions of the generators in 
A via standard functional calculus for self-adjoint operators and/or pseudo-differential operator 
calculus. 



2.3 Integrals of motion 



As the next step we shall introduce our main ansatz for the set Q of integrals of motion using the 
usual scheme of the QISM. To this aim let us assemble the generators of A into the following L- 



matrix acting on 1 



Vs. 



L xxz (w) 



i e* bu f 



s 



ueC. 



(10) 



In the definition of L xxz (it) we have used the rescaled generators e s , f s , k s which are defined by 

e s = (2sin 7 )E s , f fl = (2sin 7 ) F s , k s = K s . (11) 

Occasionally we will omit the superscript xxz for the sake of brevity. The defining relations (01 and 
© of W,(s[(2,R)) are equivalent to 



R12 (u) L 13 (u + v) L 2 3 (v) = L 23 (v) Li 3 (u + v) R 12 (u) . 
(id®A)L± = Lf s Lf 2) 



where L ± arise in the decomposition 

L(u) 

and the auxiliary R-matrix is given by 

/ sinh 7r6(w + ib) 



izbu T + 



-irbu 



(12) 
(13) 



(14) 



R(u) 



sinh 7rbu 
i sin7r6 2 e~ 7Tb 



isinirb 2 e* hu 
sinh irbu 



V 



(15) 



sinh7r6(u + ib) J 
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Out of the L-matrices we may then construct the monodromy M(it), 

M(u) = (^ U \ ^\) = Ln(«) • . . . • L 2 (u) • Li(u) . (16) 

Of particular importance is the one-parameter family of operators: 

T(u) = tr (M(u)) = A(u) + D(u) . (17) 
The trace in (fTTT i is taken over the auxiliary space, which is C 2 for the models we consider. 
Lemma 1. The operators T m which appear in the expansion 

N 

T(u) = e 1 ™" ^{~e- 27Thu ) m J m , (18) 

m=0 

are positive self-adjoint and mutually commuting, [T m , T„] = 0. 

Commutativity of T m follows by the standard argument from the relation (fT2l) : the proof of their 
positivity and self-adjointness is given in Appendixlcl 

Definition 1. Let us define the set of commuting charges as Q = {To, Ti, . . . , Tn}. 

Remark 1. Self-adjointness of the operators T m , m = 0, . . . , N ensures the existence of a joint 
spectral decomposition for the family Q. 

Let us emphasize that the crucial positivity of the operators T m is a direct consequence of the fact 
that the generators E, F and K of U q (sl(2, E)) are represented by positive operators in the repre- 
sentations Vs. This makes clear why these representations are particularly well-suited for defining 
non-compact analogues of the XXZ spin chains. We will later make a similar observation in the 
lattice sinh-Gordon model. 

2.4 Fundamental R-operator and Hamiltonian 

Our next aim is to construct a local Hamiltonian which commutes with the elements of Q. We will 
adapt the approach from IFTTl to the case at hand. The main ingredient of this approach is the so- 
called fundamental R-operator corresponding to dTOb . This operator, RJ*<? (ti), acts on V Si ®V S i 
and is supposed to satisfy the commutation relations 

R 23 (u) L 13 (u + v) L xz [v) = L lz (v) L 15 (u + v) R Zi (u) . (19) 

For our purposes it will be sufficient to deal with R(u) = R** z (u) acting on V s ®V S - 

1 The general solution R J* 3 Z (u) is needed if we wish to construct an inhomogeneous spin chain, for instance the one with 
alternating spins (see, e.g., | BD|). 
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Definition 2. Let the operator R{u) be defined by the formula 

R(u) = P w b (u + s) w b (u - s) = P D u (s) , (20) 

where P is the operator which just permutes the two tensor factors in V s ® V s , and s is the unique 
positive self-adjoint operator such that 

4cosh 2 7r6s = (n s ®7r s )A(C) (21) 

The special functions Wb(x) and D a (x) are defined in Appendix\A\ 

Theorem 1. The operator R(u) satisfies the equation ( 1791 ) where L(u) is given by MOh 

The proof of this theorem is given in Appendix iDl where the construction of the operator RJ^f (u) 
is presented for the general case, s x ^ s 2 , see equation dD.141 >. 

The operator R(u) has the following further properties 

regularity R(0) = P , (22) 
reflection property R(-u) = PR _1 (m)P, (23) 
unitarity R*(u) = R _1 (u) for u e R , (24) 

which follow from the properties of Wb(u) and D a (x) listed in Appendix lAl 

The regularity condition ( f22l allows us to apply the standard recipe [FTT| of the QISM in order to 
construct a Hamiltonian with local (nearest neighbour) interaction of sites: 

Q 

^-tr a (R n(m) • • ■ ■ • Ra2(u) ■ Ral(u)) 

N N (25) 

E—rd u D u (s n ^ n+1 ) = \ H^+i ■ 



irb 



u=0 

n—l 



We are using the following notation: We identify sn,n+i = Sn,i. the cyclic shift operator U is defined 
by U f(x X) %z, ■ ■ ■ , ^n) = f(x z , ■ ■ ■ , xn, and the subscript a stands for an auxiliary copy of the 
space V s - The trace operation is defined for an operator : V s i— > V s in the usual way: if the integral 
kernel of in the momentum representation is given by 0(k\k'), then trO = dkO(k\k). 
According to this definition we have tr a P a {, = l^. 

Substituting the integral representation (IA.20b for D u (x) into ( f25l l. we obtain the following local 
Hamiltonian density 

1 r cos(2bis n , w+ i) 
7r J sinht sinh b 2 t 

R+iO 



^f„ z +i = -- / A " • TVo, - (26) 



It may then be shown in the usual manner [FTT| that H commutes with the trace of the monodromy 
matrix, T(u), which means that 

[H, T fc ] = 0, for fc = 0,...,N. (27) 
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As in any quantum mechanical system, the fundamental problem that we would like to solve is 
the problem to determine the spectrum of H. However, thanks to the commutativity $T7\ it seems 
promising to first solve the following 

Auxiliary Spectral Problem: Find the spectrum of the operator T(u), i.e., the joint spectral de- 
composition for the family of operators Q = {To, . . . , Tn}- 

Simple counting of the degrees of freedom suggests that the spectrum of Q may be simple, i.e., that 
an eigenstate ^f t of T(u), 

T(u)*t = f(tt)*t, 

is uniquely characterized by the eigenvalue t(u) = e 7rfcNu ^m=o( — e~ 27r6u ) m t m . This would 
imply that H = H(Q), so that the solution to the Auxiliary Spectral Problem also yields the spectral 
decomposition of H. 

2.5 Classical limit 

Let us discuss the classical limit of the quantum Hamiltonian (f26b . So far we have been working in 
the units where the Planck constant H was chosen to be unity. In order to recover it explicitly, we 
have to make the following rescaling 

b —* hb , p^ft 2 p, x— »S >x, s — > 7i 2 s, s —> h 2 s, (28) 

so that we have q — e lh ' f , 7 = irb 2 . The operators e s , f s , k s are not affected by the procedure t2Si . In 
the limit h — * they become classical variables e, f , k with the following Poisson brackets obtained 
by the correspondence principle, [ , ] — > — iH{ , } 

{e,k}= 7 ke, {f,k} = - 7 kf, 2 7 {e,f} = k 2 -k- 2 . (29) 
Using the asymptotics (computed by means of contour integration) 



e -itz 

R->6 2 ./jn_jo sinh t sinh fit 2 J^ +i0 t sinh t 



h f e~ itz 1 
lim - / dt — — — = - / dt 



^ ^ T7 



(30) 



n>l 



we obtain from d26| i the corresponding classical lattice Hamiltonian density, 

H™ 1 = limMf* x „ z +1 = ilog(4cosh 2 7r6s^ n+1 ) = ±log(q 

= i log((e„f„ + i + f„e„+i) k^ x k n+ i (31) 
+ 2k r ; 2 k 2 +1 + 2cosh(27r6 S )(k r ; 2 + k 2 +1 ) + 2) . 

Here C^ 1 Tl+1 is the classical limit of the tensor Casimir operator given by (fJTJ). 
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2.6 Comparison with similar models 

L-matrix ( TTOb and R-matrix (fT3T > are suitable for the usual XXZ model as well. The only (but 
essential) difference is that in the latter case matrix coefficients of the L-matrix act on a highest 
weight module of U q (s b ) ■ We also remark that ( fTOb differs from the most commonly used "standard" 
L-matrix in that it contains extra factors e ±7rbu in the off-diagonal elements. The "standard" L- 
matrix does not satisfy (fT4l i but is symmetric (if the matrix transposition T is combined with the 
operator transposition * such that f* = e s and k* = k s ) and corresponds to the symmetric auxiliary 
R-matrix d34i >. 

Let ipb(x) denote the logarithmic derivative of the function 5&(x) defined by ( lA.15b . Properties 
( IA. 16b — (IA. 1 8b show that Sb(x) can be regarded as a b-analogue of the gamma function. The Hamil- 
tonian density d26l l rewritten in terms of ipb(x) looks as follows 



= M% + is) = ± M% - is) = ^ {M% + is) + M% - * s )) ■ (32) 



Equivalence of these expressions is due to ( IA.18I ). The last of them resembles the form of the 
Hamiltonian density of the non-compact XXX magnet [DKM| expressed in terms of the ordinary 
^-function. 

The special function Wb{u) is closely related (cf. Eq. ( ID. 71 )) to the non-compact quantum diloga- 
rithm gb{u). Counterparts of d2Qb and (ID. 14b for the compact XXZ magnet look similar in terms of 
the q-gamma function which, in turn, is closely related to the compact analogue of gb{u) given by 



It is also worth noticing that the R-operator ( f20b resembles the fundamental R-operator r (s, A) found 
in BFVII for a simpler L-matrix related to the Volterra model IIV1L The main difference is that the 
operator argument s of r(s, A) has a much simpler structure in terms of the variables p and x. It 
would be interesting to clarify the connection between these two R-operators. 



3. Lattice sinh-Gordon model 
3.1 Definition of the model 

In this section we will begin to develop the QISM for a lattice version of the sinh-Gordon model, 
which has U q (sl(2,M.)) as a quantum symmetry. We are going to keep much of the set-up from 
Section|2j but we will now be using the following L-matrix acting on C 2 <g> V a 



^(t)-aT =0 (i+^ 2ri+1 )- 




(33) 
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L-matrix d33l satisfies the intertwining relation (fT2l i where the auxiliary R-matrix is now given by 



(sinh7T&(ii + ib) \ 

sinh7r6u isin.7 | 

i sin 7 sinh irbu I 

sinh7rfe(w + i6) / 

This R-matrix possesses the following symmetry 

[R(u),a a ®a a ] = 0, (35) 

where er a , a = 1,2, 3, are the Pauli matrices. We may then proceed the same way as in Section|2]to 
define the operator T(u) and the family Q = {To, . . . , Tn} of commuting observables. We again 
find (see Appendix [CT> that the corresponding operators T m , m = 0, . . . , N are positive self-adjoint 
as a direct consequence of the positivity of e s , f s , k s . The existence of a joint spectral decomposition 
for the family Q is thereby ensured by the spectral theorem (cf. RemarkHJ. 



3.2 Fundamental R-operator and Hamiltonian 

Now our aim is to find the fundamental R-operator Rg° Sl (u) corresponding to L-matrix (1331 . For- 
tunately, it turns out that it can be constructed from the R-operator of the XXZ chain in a simple 
way. To demonstrate this, we first introduce an automorphism 9 such that 

0(p) = -p, 0(x)=-x. (36) 

It is useful to notice that 6 can be realzied as an inner automorphism, 9(0) = fiO n — 1 , where n is 
the parity operator whose action in the momentum representation is defined by (nf)(k) = f(—k). 
Notice that n is unitary and satisfies n^ 1 = si. Observe that for the representation V s we have 
(cf. ©) 

9(e) =f, 0(f) = e, 6»(k) = k 1 . (37) 
Definition 3. Let the operator Rg G s (u) be defined by the formula 

R B ° s ^(u) = (k®k)- i ^-(n®l)-Rjf s z i (M)-(l®n)-(k®k)- i ^, (38) 
where R**,? (u) is given by ( 1201 ) if s x = s z and by ( ID. 741 ) otherwise. 

Proposition 1. The operator R^ G Sl (u) satisfies the equation ( 1791 ) where L(u) is given by ( liil ). 

Proof. It will be convenient to consider R^(u) = P 21 R^ Si (u) : V Sz ® V S:L — > V Sl ® V Sl and 
analogously defined R** z (u) = P^Rff^ (u). Equation jl9[ for R|° s (u) is then equivalent to 

R SG (u) (u + v) L SG (v) = L s ° (v) Lfl (« + «) R£ («) ■ (39) 
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Let us also note that, by using (k<£> k) _1 RJ^ z (u)(k(g) k) = R xxz (u), we may rewrite the expression 
for R^°(u) which follows from d38l l as 

ft£(u) = (id®0)(ft' ia (u)), 

R' lz (u) = (k®l)~** ■ R^ z (u) • (1 k)** . ' ' 

The key to the proof of the Proposition will then be the following relation between the L-matrices 

L SG (u) and L xxz (u) 

L SG (u) = -ie~ wbs <7i k- 1 ^ L xxz {u) k l t ee CTl L'(u) . (41) 
Inserting ( flTT i into d39l l, we get an expression which contains eri L^ 3 (u) u\. Observe that (cf. d37l >) 

<7i L'(u) ai = (id <g> 0) £'(tt) . (42) 
Therefore, by using ^(OiOi) = 0(Oi)0(O 2 ), one finds that ( f39t is equivalent to 

R' 23 (u) L' 1} (u + v) L' xz {v) = L' 1} (v) L' lz (u + v) R' Z} (u) , (43) 
which is now easily reduced to the Theorem[T](its general case for R xx s z (it)). □ 

It is easy to see that properties (F22b-(l24li of the R-operator of the modular XXZ magnet hold for 
R(u) ee R^ (u) as well. Therefore, using the regularity of R(u), we can construct a Hamiltonian 
with nearest neighbour interaction of sites by using the same recipe that we used to derive ( T25l l. This 
yields 

Hl G n+ i = %d u R s n % +1 (u) # ^d®e)H™: +1 + ±\og{KK +1 ) 

u—0 ' 

1 f cos(26ts„, n+ i) v ( 44 ) 



7T 

B+iO 



7 ^ sinhtsinh^ + - l0S(k " k " +l) 



where s = (1 <g> n) s (1 (g> n) is the unique positive self-adjoint operator on "P s <8> such that 

4cosh 2 7r6s = (tt s <g> (6>o7r s ))A(C*) . (45) 

3.3 Relation with the continuum theory 

To begin with, we may first compute the classical limit, h — > 0, of d44l in the same way as we 
derived ( f3TT >. Using d37l ). we obtain 

ff^„+i = ^ log(e„e n+ i + f„f n+1 + 2(k- 1 k" 1 _ 1 + k„k„+i) 
+ 2cosh(2 7 r6s)(k j ; 1 k„ + i + k^k^)) . 

In order to establish the relation with the sinh-Gordon model it will be convenient to change vari- 
ables as follows: 

2tt6p„ = -/?$„, 4^x„ =/3(in„-$„), /3 = 678^. (47) 
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The variables $„ and Tl n will then turn out to correspond to the (discretized) sinh-Gordon field and 
its conjugate momentum, respectively. The classical field and momentum variables defined in d47b 
satisfy the Poisson-bracket relations {n„, $ m } = S nm . The Hamiltonian d46b now looks as follows 

Hn%+i = 7j log | (§ cosh f (n„ + n„ +1 ) + x -±f cosh f ($„ - $„ +1 ) (48) 

+ f cosh/3($„ - i(n„ + n„ +1 )) + f cosh /3($ n+1 - i(n„ + n n+1 )) 

+ /i coshf ($ n + $„+i) + ^cosh/3($„ + $ n+ i - i(n„ + n„ +1 ))) , 

where /i — e~ 2vbs . In order to define the relevant limit leading to the continuous sinh-Gordon 
model, let us combine the limit of vanishing lattice spacing N — * oo, A — * (i? = NA/2-7T is kept 
fixed) with the limit where the representation parameter s goes to infinity in such a way that the mass 
parameter m defined via 

imA = . (49) 

stays finite. In addition we shall assume the standard correspondence between lattice and continuous 
variables: 

n„^n(.x)A, $(a:), x = nA. (50) 

We then find the following limiting expression for the Hamiltonian density: 

Y] 3#n G „+i - const + / dx (in 2 + \{d x <$>) 2 + 2j£ cosh/3$) (51) 

thus recovering the continuous sinh-Gordon model. 

It is also instructive to see what happens to the L-matrix in this limit. In the classical continuous 
limit, i.e., when m in (|49l is kept fixed and H, A — > 0, Eqs. (|28| ) and ( TSOb show that L-matrix ( f33b 
becomes 

L SG (i)^(j;)+AC/»+0(A 2 ), (52) 

where U SG (u) is the well-known [/-matrix from the Lax pair for the classical continuous sinh- 
Gordon model ||KBl][m, 

U°°(u)-( ^ {X) isinh( W -f$(x)) \ 

Remark 2. The classical lattice Hamiltonian density d46l > resembles that found in lITall for the lattice 
sine-Gordon model. However, relation between the quantum Hamiltonians is less clear because 
the fundamental R-operator proposed in [Ta| is represented as a product of R-operators of the type 
r(s, A) which we mentioned at the end of Subsection 12.61 Possibly, recent results on factorization 
of R-operators [DKK| will help to clarify the connection between our construction and that used 
in ilTallFVirVTl . 
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3.4 Failure of the algebraic Bethe ansatz 

For the sake of clarity it may be worthwhile explaining in some detail why the algebraic Bethe ansatz 
is not suitable for the solution of the lattice sinh-Gordon model. 

To begin with, let us observe that the L-matrix d33l has no pseudo-vacuum state, i.e., a vector '5 
such that * = 0. Indeed, this would require that k„* = = k^ 1 * for all n = 1, . . . ,N. 

Such a vector does not exist. 

For the sine-Gordon model, one circumvents this difficulty by considering the composite L-matrix, 
C(u), which is product of two L-matrices IIFSTIIlKl . For 7 = 7r^, n, m e N, exponential operators 
e 1 ' 3 *™, e 4,3n " admit finite-dimensional representations. In this case there exists a vector '5 that 
is annihilated by £21 (u). This makes it possible to apply the algebraic Bethe ansatz technique. 
Let us therefore consider the analogous construction for the sinh-Gordon model. Let C SG (u) — 
L2 G (u) L\ G {u + vo), where we introduced the shift by the constant e Kin order to increase the 
generality of our consideration. We then have 

C s 2 G (u) = ie nbu (k ® e + e 7r6aJ f <g> k) - ^-^(k -1 <g> e + e _7r6ro f ® k" 1 ) . (54) 

The requirement that a vector is annihilated by £|i ( u ) is equivalent to the two equations 

(k ±x (g) e + e ±7rbro f <g> k ±1 )* = . 

We claim that there does not exist a reasonable (even in the distributional sense) state with such 
properties. Indeed, note that 

k ® e + e 71 ^ f ® k = (k^n ® 1) • (k _1 ® e + e ® k) • (k'Tn ® l)" 1 , 

where Q is the parity operation introduced in Subsection 13.21 The operator (k'^n g 1) is unitary, 
which allows us to conclude that k ® e + e^ 6 ™ f (g> k and k _1 ®e + e®k have the same spectrum. 
However, the latter operator represents (2sm r y)A(E) on V s ® V s (cf. eq. d3j). Unitarity of the 
Clebsch-Gordan maps (see Appendix ID.lt implies that this operator has the same spectrum as e s . 
The unitary transformation used in the proof of Lemma [4] in Appendix B maps e s to e 27rbx . It is 
now clear that all these operators do not have an eigenfunction with eigenvalue zero, as would be 
necessary to construct a Bethe vacuum. 

Remark 3. Keeping in mind that the sinh-Gordon variables are just linear combinations of p and x, 
cf. (|47| |. we now see quite clearly that the failure of the Bethe ansatz is connected with the fact that 
the target space (the space in which the fields take their values) is non-compact. We expect this to 
be a general lesson. 

4. Q-operator and Baxter equation 

As an important first step towards the solution of the Auxiliary Spectral Problem we shall now find 
necessary conditions for a function t(u) to be eigenvalue of the operator T(it). In order to do this 
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we are going to construct an operator Q(u) which satisfies the following properties: 

(i) Q(u) is a normal operator, Q(u)Q* (v) = Q* (v)Q(u) , 

(ii) Q(u)Q(») = Q(w)Q(«), 

(iii) Q(u)T(«) = T(u)Q(«), (55) 

(iv) Q(«) T(«) = (a(u)) N Q(« - ib) + (d(u)) N Q(u + ib). 

The first and the second property imply that all operators Q(u), u E C can be simultaneously 
diagonalized and their eigenvectors form a complete system of states in the Hilbert space. The third 
and the fourth property imply that T(u) will be diagonal whenever Q(u) is. One may therefore 
consider the spectral problem for Q(u) as a refinement of the spectral problem for T(u). 

Let us now consider an eigenstate *f? t for T(u) with eigenvalue t(u), T(u)\I/ t = t(u)^f Thanks to 
property (iii) above we may assume that it is an eigenstate for Q(u) as well, 

Q(u)*t = *(«)**■ (56) 
It follows from property (iv) that the eigenvalue q t (u) must satisfy the so-called Baxter equation 

(57) 



t(u)q t {u) = (a(w)) N q t (u ~ ib) + (d(u)) q t (u + ib) . 



We will construct the operator Q(u) explicitly — see Subsection |4.1| This will allow us to determine 
the analytic and - for the lattice sinh-Gordon model with N odd (the SGo-model) - the asymptotic 
properties that the eigenvalues qt(u) must have, namely 



(i) qt(u) is meromorphic in C, with poles of maximal order N in T_ S UT S , 
where T s = {s + i(f +n6 + m6" 1 ), n,meZ- }, T s = (T s )* , 
exp ( + niN(s + |Q)u) for |u| — > oo, |arg(u)| < |, 
exp ( — 7T?N(s + |Q)m) for \u\ — > oo, |arg(u)| > |. 



(ii) <r» 



(58) 



The derivation of these properties is discussed in Subsection 14.21 This means that there is the fol- 
lowing necessary condition for a polynomial t(u) to be eigenvalue of the operator T(u): t(u) can 
only be an eigenvalue of T(u) if there exists a meromorphic function qt(u) with singular behavior 
and asymptotic behavior given in ( 1581 which is related to t(u) by the Baxter equation ( f57b . 

The problem to classify the solutions to this condition is of course still rather nontrivial. However, 
previous experience from other integrable models suggests that the Baxter equation supplemented 
by the analytic and asymptotic properties ( 1581 is indeed a useful starting point for the determination 
of the spectrum of the model, see also our Subsections 16.11 and 16.21 for some further remarks. We 
will discuss in the next section how the separation of variables method may allow us to show that 
the conditions above are also sufficient for t(u) to be an eigenvalue of T(it). 

Convention: We will use the superscripts sg and xxz to distinguish analogous operators within the 
two models we consider. However, we will simply omit these superscripts in any equation which 
holds in the two cases alike. 
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4.1 Explicit form of Q(u) 

Let us now describe explicitly the Q-operators for the models that we introduced in Sections [2] 
and[3] For this purpose we will work in the representation where the operators x r , r = 1, . . . , N are 
diagonal. This representation will be called the Schrddinger representation for the Hilbert space of 
a lattice model. Let x = (x x , . . . , x N ), x' = (x' x , ... ,x' n ). We will denote the integral kernel of 
the operator Q(u) in the Schrddinger representation by Q„(x, x'). We will also use the following 
notations 

a = s + |Q , a = s - \Q , (59) 
where s stands for the spin of the representation V s - 

Definition 4. Let the Q-operators (u), b = xxz, sg, be defined in the Schrddinger representation 
by the following kernels 



Qt itt (x,x') = (60) 

N 

= (D_ s (u)) J| Di (s _ u) (x r - x' r )Di^ +u) (x r -l ~ £], x' r ) D- S (x r - S b X r -l) , 



r=l 
N 



Qi ;u (x, x') = J| Di M (i r ~ O £>_ i („ +0 -)(»r " £b A,(a> ~ £b , (61) 

r=l 

where s xxz = 1, e SG = — 1, ant/, in f/ie sinh-Gordon case, s is related to the parameters m, A as 
in dH. 

Theorem 2. Lef T (u), b = xxz, sg be the transfer— matrices corresponding to the L— matrices ( 1701 ) 
one/ ( I35I ). 

(i) The operators (u) satisfy all relations in ( |55| ). 

f iij 77ze Baxter equation ( 1551 -iv) holds for Q j_ (u) wif/i the following coefficients 

a xxz (u) = 2 sinh 7rfe(u - cr) , d xxz (w) = 2 sinh tt&(w + a) , 
d SG (u) = a SG (-u) = e 7rb(u+l ^ + {m^f e -^+m . ,6 ~ ) 

(iii) 77ie operators (u) satisfy the relation 

Q h +{u)Qt_(v) = Ql{v)Qt + (u). (63) 

The proof of this Theorem is given in Appendix [E] It is worth noting that Q_(u) and Q+(u) are 
related by hermitian conjugation as follows. 

Q-(u) = (-D_s(?i)) N (Q+(u))* . (64) 
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This allows us to mostly focus on Q b (u) = Ql(u), but it is nevertheless sometimes useful to 
consider Q^.(u) as well. The corresponding eigenvalues qf{u) are consequently related as 



<£-(«) = (£-(«)) q+(u). (65) 



Relation ( I65l l will imply that and q t have the same analytic and asymptotic properties i 

Remark 4. It is sometimes useful to observe (see Subsection \E.2\ that the operator Q b (u) can be 
factorized as follows: 

Q\u)=Y b (u)-Z, (66) 
The operators Y b (u) and Z in d66*l l are represented by the kernels 

N 

^«( x > x ') = Yl D ±(u-*)( x r - e b x' r+1 ) D_i {u+a) {x r - x' r ) , (67) 

r=l 

N 

Z(x,x') - (io 6 (*f - 2s)) N J] D s (x r - x' r ) , (68) 

r=l 

respectively, where e X xz = 1, £sg = — 1- 

Remark 5. The relations (f55t do not define the Q-operator Q(u) uniquely. For instance, for a given 
Q(u), relations d55l > are also fulfilled for Q'(u) — (<£>(u)) Q(u), where ip(u) is a scalar function 
and is a unitary operator that commutes with Q(u) and T(u). The coefficients in Baxter equation 
(fSBT-iv) for Q(u) and Q'(u) are related via 

a'(u) = o(«) , = d(u) . (69) 

ip\u — ib) (p[u + lb) 

Thus, there is no canonical way to fix these coefficients. However, their combination a(u)d(u — ib) 
remains invariant; its value is related to the quantum determinant if the latter can be defined for the 
L-matrix of the model in question (see Appendix lC.2l ). 



4.2 Analytic properties of eigenvalues of Q(u) 

We now turn to the derivation of the analytic properties of eigenvalues of Q(u). More precisely we 
shall prove the following: 

Theorem 3. 

( ;) The operators Q^. xz (u) and Q j_ G (u) are me romorphic functions ofu in C with poles of maximal 
order N contained in T_ s U T s . 

(ii) Denote Q j_ Go (u) = Q s ± (u) for N odd. This operator has the following asymptotic behavior 

/~vSGo / \ f Q+oc cxp( + i7rNcru) for |u| -> oo, |arg(u)| < f , 
Q + (u) ~ < (70) 
I Q_ so exp ( — iirNau) for |tt| — > oo, |arg(u)|> : |, 



20 



where Q ±00 are commuting unitary operators related to each other as follows 

Q_=QQ +oc . (71) 
Here Q is the parity operator ( its action in the Schrddinger representation is given by ( l<S2| ) j. 
(Hi) Ql Go (u) has asymptotic behavior of the same form ( 1 701 ) with Q ioo replaced by Q*^,. 

Proof. In order to prove part (i) of the Theorem it clearly suffices to consider the corresponding 
statements for the eigenvalues qf(u). Let us first explain why the properties of qf{u) and q^(u) 
described in the theorem are the same. Given that the theorem holds for (u), we infer that poles of 
qt(u) are contained in T S UT_ S . But, since they are of maximal order N, they cancel in d65l ) against 
the N-th order zeroes of (-D_ s (u)) N (see properties of D a (x) in Appendix IA. 31 ). Thus, the only 
possible poles of <j£~ (u) are those of (D_ s (u)) N , i.e., they are of maximal order N and contained in 
T_ s U T s . 

The proof of part (i) will be exactly analogous for the cases of the XXZ magnet and the sinh-Gordon 
model. Only the latter case will therefore be discussed explicitly. We will study the equation ( 1561 ), 
which is equivalent to 

/ dxdx'$(x)Q^ u (x,x')* t (x') = g+(u)<$|*,>, (72) 

JR2N 

for some test-function ^(x) £ 7^® N , where % is the space of test-functions which is canonically 
associated to the representations V s as shown in Appendix|B] In order to find the analytic properties 
of qf (u) let us use Q s +(u) = Y SG (u)Z to represent the left hand side of d72| > as ( <E>' | Z^ t ), where 

$V) = / rfx$(x)K u SG (x,x'). (73) 

With the help of the Paley- Wiener theorems one easily finds that the condition ^(x) G implies 
that ^(x) is entire analytic w.r.t. each variable Xk and decays exponentially as 

|$(x)| - e -*Q\ x >*\ for \x k \ -> oo. (74) 

The kernel y,f G (x, x') has the same asymptotics w.r.t. its x k variables, as seen from equation j67l ) 
and relation JA.25b . Therefore the convergence of the integrals does not represent any problem. 
Combined with the observation that the left hand side of d73l is the convolution of two meromorphic 
functions we conclude that the only source of singular behavior is the possibility that the contours 
of integration in (l73b may become pinched between poles of the integrand approaching the contour 
from the upper and lower half planes, respectively. With the help of dA. lib one easily compiles a list 
of the relevant poles of the kernel Y£ a (x, x') as given in (1671 : 

Upper half plane H_ : (1) x r € x' r — j(u + a) + T , 

(2) x r € -x' r+1 + \{u - a) + T , 

Lower half plane H + : (l') x r 6 x' r + ^(u + a) — To, 

(2') x r G -x' r+1 -\{u-o)- T . 
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Pinching of the contour between poles from the upper and lower half planes would produce the 
following series of poles: 

(11') u + se+To, (12') 4 + x' r+1 - s G -T , 

(22') u-se-To, (21') 4 + x' r+1 + s € +T . ! 

We observe in particular that none of the poles of $'(x') happens to lie on the real axis, which rep- 
resents the contour of integration for each of the integrals over the variables x' k in (1721 1, Taking into 
account the exponential decay of Yjf G (x, x') for \x' k \ — » oo, we may conclude that the integration 
over x' in ( 1721 converges nicely. It follows that the left hand side of (l72l defines a meromorphic 
function of u with poles listed on the l.h.s. of (175) . 

In order to verify part (ii) of the theorem let us note that Yjf G (x, x') has the asymptotic behavior 

K u SG (x,x') ~ TT e T 2 ^(^ +1 +<) for |u| - oo, ( |arg(u)l < f ' (76) 

r=i I Iarg(u)| > f • 

as follows straightforwardly from dA.26) . In order to check that the integral obtained by exchanging 
the limit for u — > ±oo with the integrations in d72l ) is convergent let us note that performing the 
integration over the variables x r yields the Fourier transformation "^(k) of $(x), with argument 
k = (x^ + x' z , . . . ,x' N + x'-J. Note that the change of variables k = k(x') is invertible for N odd. 
We may therefore represent the integration over x' by an integration over k. The nice asymptotic 
properties of 4>(k) which follow from our requirement $ € (T s ) ® N ensure the convergence of the 
resulting integrals. Part (ii) of the theorem therefore follows from (|76| | and d66l l. 

The proof of part (iii) of the theorem is immediate, if d64l i and (IA.25) are taken into account. □ 

Remark 6. Let us comment on the nature of the problems which prevented us to determine the 
asymptotics of the Q-operators in the remaining cases. In both remaining cases one must observe 
that the change of variables k = k(x') is not invertible, which implies divergence of the integral 
over x'. This is closely related to the fact that the leading asymptotics of T(u) for \u\ — * oo 
introduces a quasi-momentum p which has purely continuous spectrum (see equations (IC.13b . and 
( IC. 1 5b in Appendix0. It follows that T(u) can not have any normalizable eigenstate. Instead one 
should work with the spectral representation for p , H ~ L dp TC Po where the elements of H Vo 
are represented by wave-functions of the form fy Po (x) = e 2mp ° x ° $(in-i — Xo, . ■ ■ ,Xi— x ). This 
seems to complicate the analysis considerably. We nevertheless expect results similar to (|58l-ii) to 
hold for the remaining cases as well. 

4.3 Self-duality and quantum Wronskian relation 

The explicit form (l60ll-(l6TTi of the Q-operators along with the properties of D a (x) listed in Ap- 
pendix [A3] show that Q±(u) are self-dual with respect to the replacement b — > b^ 1 . Therefore, the 
Q-operators also satisfy the dual Baxter equations, 

J\u) ■ Q b ± {u) = (a'{u)f Q b ± {u - lb' 1 ) + (dV)) N Q±(« + lb' 1 ) , (77) 
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where T b (u), b = xxz, sg, denote the transfer-matrices corresponding to the L-matrices ( [Tol l and 
d33l i with b replaced by b^ 1 . These are the transfer-matrices of the modular XXZ magnet and lattice 
sinh-Gordon model with Uq(sl(2, R)) symmetry, where q — e Z7rb . The coefficients a(u), d(u) in 
(TTTb are similarly obtained from those in (l62b by the replacement b — ► In the sinh-Gordon case 
the mass m b -i is related to the representation parameter s via \m\ /&A = e _7rf ' s . 

This self-duality has remarkable consequences, which we shall work out explicitly for the case of 
the sinh-Gordon model with odd N. We will take advantage of the freedom pointed out in Remark|5] 
to renormalize the operator Q SG (u) as follows: 

Q(u) = Q* +0O e-^ N ( u2+a2+5 +)Q s + G {u), (78) 

where 5 + = |(6 + b^ 1 ) and Q+oo is the unitary operator which appears in the asymptotics ( fTOl i. 
The Baxter equation for Q(it) will then take the following form: 

T SG (u) • Q(u) = (a{u)f Q(u - ib) + {d{ujf Q(u + ib) , 
where d(u) = a(-u) = 1 + (^) 2 e "^K2«+») . 

The normalization of the operator Q(u) was for later convenience chosen in such a way that Q(u) ~ 

e±nm -,fN(„ 2 + a 2 + ij) . j forRe ( u ) ^ ±ODi 

Theorem 4. The operator Q(u) fulfills the following quantum Wronskian relation: 

Q(u + iS+)Q{u - i6+) - Q(u + iS-)Q(u - iS-) = W N {u)-l. (80) 
where W N {u) = e -"N(u 2 + ff 2 ) ( Da ( u f) ~ N and S ± = ± b). 

Proof. Let W(m) be the left hand side of d80l . A straightforward calculation, using the Baxter 
equation and its dual form, shows that \N(u) satisfies the following two functional relations: 

v 2 V coshirb^iu + a) J K 2 ; 

A solution to both functional relations is given by the expression on the right hand side of d80b . If 
b is irrational it suffices to notice that W(u) is meromorphic in order to conclude that the solution 
to the system dBTT i must be unique up to multiplication by an operator which does not depend on u. 
This freedom can be fixed by comparing the asymptotics of both sides for Re(it) — > ±oo using 
equations d70b.(F7TT> (for considering the asymptotics as Rc(u) — * — oo, it is helpful to notice that 
Q 2 = 1). In order to cover the case of rational b let us notice that both sides of the relation ( f80b can 
be analytically continued from irrational values of b to the case where b is rational. □ 

From the proof of this theorem it is clear that the main ingredient is the self-duality of our repre- 
sentations V s . We therefore expect that a similar result will hold for the remaining cases as well. 
However, at present we do not control the asymptotics of the Q-operators sufficiently well in these 
cases. 
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4.4 Parity and cyclic shift 

Let us consider the cyclic shift and parity operators U and ft defined respectively by 

U f(x x ,x z , . . . , en) = f(x z , ...,xtx,x x ), ft f(x x , . . . , x N ) = f(-x x , -x N ) . (82) 

These operators commute with each other and also with Q±(w) (as can be easily seen from d60l>— 
doTh. Hence they must commute with T±(u), as can also be verified directl)Q It follows that 
eigenstates ty t may be assumed to be simultaneously eigenstates of U and ft, 

ft* t = ±* t , U* t = e 2 "t^ t7 m=l,...,N. (83) 

It is therefore useful to observe that U and ft can be recovered from Q±(u) as follows. First one 
may notice that the integral kernel d60b simplifies for the special values u = ±a thanks to the 
relation ( 1A.30I ). Explicitely, we have 

Q$» = (u7 6 (2s+iQ))" N -l, (?+(-*) = (w 6 (2 S +fg))~ N -U- 1 -ft b , (84) 

where ft b is defined in dE.321 >. Combining this observation with d83l . we conclude that 

q+(a) = (w b (2s+±Q)y N , q+(-a) = ±e 2 ^ (w^s+^Q))^ , (85) 

where the — sign in the second expression can occur only in the sinh-Gordon model. 

Secondly let us observe that eqs. d65l > and d85l l imply that u = ±<r are poles of order N for q^(u). 
Indeed, using doTi ( IA.30I ), and ( IA.12I I. we find 

Q b _(a + e) ~ (^) N • 1 , + e) ~ (^) * ■ ft b ■ U , (86) 

as e — > 0. This implies 

?r(* + <0 - (^) N , ?r(-* + <0 ~ ±e 2 "- (^) N , (87) 

where the — sign in the second expression can occur again only in the sinh-Gordon model. 

5. Separation of variables 

In the previous section we have identified necessary conditions for a function t(u) to be an eigenvalue 
of T(u). If we were able to show that these conditions are also sufficient, we would have arrived at 
a useful reformulation of the Auxiliary Spectral Problem. 

A promising approach to this problem is offered by the separation of variables method pioneered 
by Sklyanin BS2I IS31 . The basic idea is to introduce a representation for the Hilbert space H of 

2 Indeed, for U this is obvious from the definition (TTJ and for ft it follows from the observation that ft L xxz (u) ft = 
cri e - nburT 3 L xxz («) e wbua 3 oi and ft L SG (u) ft = a x L SG (u) cti . 
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the model in which the off-diagonal element of the monodromy matrix M(u), the operator B(u), is 
diagonal. 

For simplicity of exposition let us temporarily restrict attention to the case of the sinh-Gordon model 
with N odd. The operator B(u) has the following form: 

N 

B(«) = -i N e N7rfc («-s) (~) m e - 2m7rbu B m . (88) 

m=0 

By Lemma|5] the operators B m , m — 0, . . . , N are positive self-adjoint. Basic for the separation of 
variables method is the validity of the following conjecture. 

Conjecture 1. The joint spectrum of the family of operators {B m ; m = 0, . . . , N} is simple. This 
means that eigenstates of B(u) are uniquely parameterized by the corresponding eigenvalue b{u). 

This conjecture can be supported by counting the degrees of freedom. However, it is not easy to 
provide a rigorous proof (see also Remark[7]below). 

The function e~ wbNv, b(u) is a polynomial in the variable A = e - 27rbu . It can conveniently be repre- 
sented in the following form 

N 

b(u) = -(2z) N e-^ 8 Y[ smhirb(u - y k ) . (89) 

fc=i 

The variables y k , k = 1, . . . , N are uniquely defined up to permutations once we adopt the conven- 
tion that lm(y k ) G (—55, 55]- This means that the representation for the Hilbert space H in which 
B(it) is diagonal may be described by wave-functions ^(y), y = (y 1 , ■ ■ ■ , y N ). This representation 
for the vectors in H will subsequently be referred to as the SOV representation. 

We will then show that the Auxiliary Spectral Problem, T(u)^!t — f (w)^tj g ets transformed into 
the system of Baxter equations 

t(y fc )*(y)= [(a(y k )) N T^ + (d(y k )fT+]^(y), fe=l,...,N, (90) 

where the operators Tfe are shift operators defined as 

T±*(y) = *( yi ,...,y k ±ib,..., m ). (91) 

The coefficients in front of the shift operators in d90T l depend only on a single variable y k , which is 
the crucial simplification that is gained by working in the representation where B(u) is diagonal. 

The key observation to be made at this point is that the same finite difference equation (1901 was 
found in the previous Section |4]in connection with the necessary conditions for a function t(u) to 
represent a point in the spectrum. It now remains to observe that any function q t (u) that fulfills the 
necessary conditions d57li,(l58ll can be used to construct 

N 

*t(y) = II (92) 

k=l 
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The fact that d92l defines an eigenstate of T(u) is verified by comparing d90l l with (T5Tb . The main 
point that needs to be verified is whether the function ty t (y) actually represents an element of Tt , 
i.e. whether it has finite norm. The scalar product of vectors in 7i can be represented in the form 

(*z|*0 = f dfi(y) (* a |y)(y|* a ) (93) 

We clearly need to know both the range Y of values y that we need to integrate over, as well as the 
measure dfi(y) of integration to be used. We will be able to determine the measure rf/i(y) provided 
that the following conjecture is true. 

Conjecture 2. The functions b(u) of the product form ( 1891 ) that describe the spectrum ofB(u) have 
only real roots, i.e., y^ G M for k = 1, . . . , N. 

We will discuss the status of this conjecture after having explained its consequences. Conjecture [2] 
directly implies that Y = R N in d93l ). Assuming the validity of Conjecture [2] we will show in 
Proposition|2]that the measure dfi(y) can be represented in the following explicit form 

N 

rf A*(y) = II d Vk Y\_ 4sinh7r6(y fe - yi) sinh7r6 _1 (j/ fe -yi). (94) 

k=l Kk 

Knowing explicitly how to represent the scalar product of fi in the SOV representation finally allows 
us to check that any solution of the necessary conditions ( 1571 ), (1581) defines an eigenvector | \l/ t ) of 
T(it) via ([92]). In other words: The conditions d571 l. d58l ) are not only necessary but also sufficient 
for t(u) to be an eigenvalue of a vector \^t) E T~t- 

Remark 7. Our claim that the conditions (l57l i. ( 1581 are also sufficient for a function t(u) to repre- 
sent a point in the spectrum of T(u) does not seem to depend very strongly on the validity of the 
Conjecture[2] In this sense the conjecture mainly serves us to simplify the exposition. 

In any case it is a problem of fundamental importance for the separation of variable method to 
determine the spectrum of B(u) precisely. Even in simpler models which have been studied along 
similar lines like the Toda chain IIKL1 or the XXX spin chains [DKM | there does not seem to exist a 
rigorous proof of the analogous statements. The explicit construction of the eigenfunctions of B(u), 
which may proceed along similar lines as followed for the Toda chain in |KL| or for the XXX chain 
in llDK M I. should provide us with the basis for a future proof of Conjecture|2]or some modification 
thereof. 

Remark 8. Within the Separation of Variables method the Auxiliary Spectral Problem gets trans- 
formed into the separated Baxter equations d90l ). However, these finite difference equations will 
generically have many solutions that do not correspond to eigenstates of T(w). 

In order to draw a useful analogy let us compare the situation with the spectral problem for a differ- 
ential operator like h = —dy + V(y). One generically has two linearly independent solutions to the 
second order differential equation like {—dy + V(y))ip — £ip for any choice of £. The spectrum of 
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h is determined by restricting attention to the subset of square-integrable solutions within the set of 
all solutions to the eigenvalue equation. 

From this point of view we may identify the conditions ( f58l > on analyticity and asymptotics of the 
function q t (u) as the quantization conditions which single out the subset which constitutes the spec- 
trum of T(u) among the set of all solutions of (l90l l. 

5.1 Operator zeros of B(u) 

The adaption of Sklyanin's observation to the case at hand is based on the following observations. 
First, by Lemma|5] the operators B m introduced in d88l are positive self-adjoint. Taking into account 
the mutual commutativity (which follows from ( 1 100b ). 

[B m ,B„] = 0, (95) 

leads us to conclude that the family of operators {B m ; m = 1, . . . N} can be simultaneously diago- 
nalizecH. 

Conjecture Q] implies that the spectral representation for the family {B m ;m = 1, . . . N} can be 
written in the form 

|*) = / d!/(b)|b)(b|*) (96) 

where | b ) is a (generalized) eigenvector of B m with eigenvalue b m , and we have assembled the 
eigenvalues into the vector h — (b 1 , . . . ,b N ). 

It now turns out to be particularly useful to parameterize the polynomial of eigenvalues b(u) in terms 
of its roots. This representation may always be written as follows 

N 

b(u) = b(u\y) = -(2i f e-* ms J] sinli7r6( U - y k ) , (97) 

fc=i 

where y = (y 1 , . . . , y N ). The variables yk are either real or they come in pairs related by complex 
conjugation. The variables are uniquely defined up to permutations if one requires that lm(yk) € 
(— Jg]. We will assume that yk E M. according to Conjecture|2] It then follows that the spectral 
representation ( t96b can be rewritten as 

|*) = f d|i(y)|y)(y|*). (98) 

However, points y, y' in M N which are obtained from each other by the permutation yk <-* yi will 
correspond to the same eigenstate of B(u). This means that the spectral representation for B(u) can 
be used to define an isomorphism 

H ~ L 2 (M N ; d/i) Symm , (99) 

3 By Lemma|6] we have Bq = (Bn) — 1 . 
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where L 2 (R N ; d/i) Symm is the subspace within L 2 (R N ;d/j,) which consists of totally symmetric 
wave-functions. 

Despite the fact that H is isomorphic only to a subspace in L 2 (R N ; d/i) it will turn out to be useful 
to extend the definition of the operators A(u), B(it), C(u), D(u) from L 2 (R N ; G?/i) Symm , where it is 
canonically defined via d99b to L 2 (R N ; dp). As a first step let us introduce the operators which 
act as yfe | y ) = ?/fe | y ). Substituting y k — > yfe in d97l l leads to a representation of the operator B[u) 
in terms of its operators zeros y k - 

5.2 Operators A(u) and D(u) 

Monodromy matrices of the modular XXZ magnet and the lattice sinh-Gordon model satisfy the 
exchange relations (fT2l . where the R-matrix is given by ( fT5l ) and d34b . respectively. Among these 
relations we have, in particular, the following 

[B(t*),B(i;)]=0 ) [A(u),A(«)]=0, [D(u),D(«)]=0, (100) 
smhirb(u-v + ib)B(u)A(v) (101) 

= swhirb(u - u) A(v) B(u) + i? 2 b 3 (u - v) B(v) A(u) , 
sinh7r6(u-u + i6)B(«)D(u) (102) 

= sinli7r6(u - v) D(u) B(v) + i? 3 b 2 (u - v) B(u) D(v) , 

^23 XZ («) - R% Z {-u) = ^ bu sin 7 , = = iS in 7 . (103) 

We are now going to show that there is an essentially unique representation of the commutation 
relations ( 1 1 00b — (fT02b on wave-functions ^(y) which is such that B(u) is represented as operator of 
multiplication by b(u) = b(u\y), cf. d97l i. 

To this aim let us consider for fc > 1 the following distributions: 

(y|A(|fe) = lim (y|A(«), (y| D(y k ) = lim (y| D(u) . 

u—>yk: u^y k 

These distributions, which will be defined on suitable dense subspaces of L 2 (R N , dp), can be re- 
garded as the result of action on (y| by operators A{y k ) and D(yfc) with operator arguments sub- 
stituted into iC3\ . ( IC.5I ), ( IC.6I ), ( IC.81 l from the left. The commutation relations (1 1 lb — ( 1 1 02b imply 
that (y| A(y k ) and (y| D(y k ) are eigenstates of B(u) with eigenvalues b(u\y'), with y' k — i/kT ib, 
respectively, y', = yi otherwise. This, along with relations ( 1 100b . leads to the conclusion that the 
action of the operators A(y k ) and D(y^) on wave-functions W(y) = (y|^) can be represented in 
the form 

A(y fc )*(y) = aN,fc(2/fe)T-*(y), D(y fc )*(y) - d N , k (y k ) T+*(y) , (104) 
where T^ are the shift operators defined in equation d9TT ). The functions a^, k (yk) and d^ tk (y k ) are 
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further restricted by the following identities: 

det 9 M SG (u) = A SG (u) D SG (u - ib) - B SG {u) C SG {u - ib) (105) 
= (4e- 27rbs cosh7r6(s + u - £§) coshvr&(s -u + i|)) N • (106) 



These identities are proven in Appendix IC.2I It follows then from ( 1104t that det 9 M(y k ) — 
ON,fe(2/fe)rfN I fe(?/fc — ib). Not having specified the measure /i(y) yet leaves us the freedom to multiply 
all wave-functions ^(y) by functions of the form J| fc fk{yk)- This allows us to choose 

aN,k(Vk) = (a(yfc)) N > d N ^ k (y k ) = (d(y k )) N , (107) 

where a SG (u) = d SG (-u) = e~* bs 2coshvr6(M - s - i|) 

no8) 

_ e -7rfe(«— i|) _|_ ^mA^ 2 e 7rb(u-i|) 

We have used that the sinh-Gordon parameters m, A are related to s as in (|49l . 

In the special case that ^(y) is an eigenfunction of the transfer-matrix T(u) with eigenvalue t(u) 
we get the Baxter equations ( |90l from T(yfc) = A(y/-) + D(y/ £ ) and equations dl04t and dl07| ), as 
advertised. 

It will be useful for us to have explicit formulae for A(u) and D(u) in terms of the operators y k 
and T^ . In the case of N odd we may use the following formulae: 

k=l l^k w y ' 

fe=i ;#fc K,K yiJ 

These formulae are easily verified by noting that the number of variables y^ coincides with the 
number of coefficients in the expansion (IC.6b and ( IC.8b - It follows that the polynomials A SG (u) and 
D SG (u) are uniquely determined by their values A SG (y k ) and D SG (y k ), k = 1, . . . , N. 

5.3 Sklyanin measure 

We furthermore know that the operators A m and D m which are defined by the expansion 

N-l 



A SG (u) = _ j e -^jN e 7r6N(«- 5 ) J2(-) m e- 2m7rbu A m , (111) 

m— 

N-l 

D SG (u) = -ie" 7 ^ e ^ bN («- s ) ^ (-) m ^ 2 " ta D m (112) 



m=0 

are positive (Lemma[5]in AppendixlO. 
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Proposition 2. There exists a unique measure dp,(y) such that the operators A m and D m on 
L 2 (R N ; d/i) are positive. This measure d[i can be represented explicitly as 

N 

d Ky) = Yi dyk II 4sinll7r& (^fe - yi) sinh7r6 _1 (y fc - y t ) . (113) 

k=l Kk 

Proof. The similarity transformation ^(y) = Xa(y) ^(y)> where 

N 

x*(y) = II (^'^(v* - s))- N JJK(2/fc - + , (114) 

k=l Kk 

maps to a representation in which the operator A SG (u) is represented as 

N 

A SG (u) = II sinh7r6(u-y z )Tfc . (115) 

fc=i ;^fc 

Expanding in powers of e^ 6 " yields a representation for the coefficients A^ 3 that appear in the 
expansion JC.3I ) which takes the form 

N 

Am = 5^Pmfc(y)T^ , (116) 

fc=l 

The coefficients p m k{y) in dl 16t are positive for all y 6 M N , and p m fe(y) does not depend on y^. 
We are next going to show that the positivity of implies that T^ must be a positive operator in 
L 2 (M N ; dp.). Let us keep in mind that T^ satisfies the commutation relations 



e 



itnj- e m = e w*«T-. (117) 



If there was any negative contribution to the expectation value ($ | J\=i Pmk(y)T fe | we could 
make it arbitrarily large by means of the unitary transformation |$) — > e ltyfc |<I>). It follows that 
<* I Pmfe(y)Tfe | *> > for any k = 1, . . . , N. 

It remains to notice that, since the p m k (y) are non-vanishing, vectors of the form \Jp m k{y) | $ ) 
form a dense subset in L 2 (R N ; dp) . This finally allows us to conclude that must be a positive op- 
erator. But this furthermore implies that (T^) 4 ' is a unitary operator which satisfies the commutation 
relations 

(j-yl e ity m = ex p(iWt4 m ) e itym (T^) 2 ' . 

It is well-known that the representation of these commutation relations by unitary operators is 
essentially unique. The measure which defines the corresponding Hilbert space is just dv{y) = 

nLi d vk- 

It remains to return in ( |98l l to the original representation via (II 141 >. that is to compute 
dfJ-(y) = IXA(y)| 2 dv{y). Using ( IA.7ft — ( tATTOb to simplify the resulting expression yields the for- 
mula for dp(y) stated in Proposition^ 
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For the operator D SG (u), a completely analogous consideration applies with the transformation 
*(y) = X D (y)$(y), where 



(y) = ft ( e ~™ ykS MVk + «)) N l[(w b (y k - Vl + |Q)) 1 . (118) 



) N 

k=l l<k 



Thanks to the relation (lA.lOt and the Conjecture|2] we have |xA(y)| 2 = IXo(y)| 2 - This leads to the 
same measure dfi(y) given by dl 131 >. □ 



5.4 Remaining cases 

To end this section let us briefly discuss the necessary modification in the cases of the modular XXZ 
magnet and the lattice sinh-Gordon model with even N. The main new feature that arises in these 
cases is the existence of a quasi-momentum y which first appears in the expansions 

N-l 



b xxz (u) = 2 N - 1 i e 7rfc ("+y°) JJ swh<irb(u - y k ) , (119) 

fe=i 

N-l 

b SG (u) = -(2i f~ 1 e^ y "' Ns '> Y[ sinhirb(u-y k ), N - even , (120) 



The variable y requires a slightly different treatment compared to the y k , k > 1. Considering the 
states (where n = for the modular magnet and k = i/2 + bs for the sinh-Gordon model with 
even N) 

(y|A ^ lim e^ K -^(y\A(u), 
(y|D ^ lim e *N(«-6u) ( y | D (u) , 

u— »-j-oo 

(y| An = lim (-fe^+ b ^{y\A(u), 

u—> — oo 

(y| Dn = Hm (_)N e -N( K+ 6„) (y| D(u) > 



and taking into account the asymptotic behaviour of B(u) and of the coefficients ( |103t at u — ► ±oo, 
we infer from the relations d 1 1 b — ( fT02l > that 

A *(y) = a (y )T+*(y), D *(y) = d {y ) Tj*(y) , 
A N *(y) = a N (y ) T~*(y) , D N *(y) = d N {y ) T+*(y) , 



where the shift operators T^ are defined analogously to (|9TT l for the variable y . It follows from 
dC13]l and (|C15jl that we have a (y) = d$(y), d (y) = a N (y), a {y)a N (y + ib) = d (y)d N (y - 
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ib) = 1. Noting that 

*(y)= lim r e ir6N(2 ' ,±(a - au)) det < ,M(u)l« f (y)= (121) 

f A D #(y) = a {y - ib)d (y )^(y), u -> +oo; 
[ A N D N *(y) = a N (y + ib)d N (y )9(y), u -> -oo , 

allows us to choose 

a = d = on = c^n = 1 ■ (122) 
The resulting equation for the y -dependence may therefore be written as 

lim [(±)N e -N( KT M t(u) ] $w = (T + + T -) $(y) _ (123) 

This relation supplements the Baxter equations d90b in the cases of the modular XXZ magnet and 
the sinh-Gordon model with even N. 

In the case of the modular XXZ magnet we furthermore find a small modification in the form of the 
coefficient functions a(u) and d(u) which appear in the Baxter equations. These follow from the 
following formula for the q-determinant (see Appendix lC.2l ): 

det 9 M xxz (u) = A(u) D(u - ib) - q' 1 B(u) C{u - ib) (124) 
= (4 cosh7r&(s + u — i|) cosh7r6(s — u + i§)) N , (125) 

The resulting expressions for a xxz (w) and g? xxz (u) will be 

a xxz (u) = d xxz (-u) = -2i coshnb{u - s - i|) . (126) 

The existence of the "zero mode" y also leads to modifications in the formulae for A(it) and D(u). 
For the modular magnet and the sinh-Gordon model with even N, the number of coefficients in the 
expansion iCM and (IC.8b exceeds by two the number of the operators y/-. However, in these cases 
we know the asymptotics of A (it) and D(u) and therefore we will need the following interpolation 
formula. 

Lemma 2. Let e N7Tbu P(u) be a polynomial in e 27Tbu such that 
P(u) 



e +7Tb(N M + Po ) ; faru^+oo. 
(_)N e -*6(N«+p.) f foru __ OQ 



For an arbitrary set of variables y\, . . . , yjq-i such that yk ^ yifor all k ^ I we may then write 

N-i N-i ' h bf - ) 

P(u) = smhTTb(u +p + p N ) TT sinh7r6(u - y k ) + V TT Sm ^ Vl> P{y k ) , 

where p N = J2k=i Vk- 



32 



Thus, for the modular magnet and the sinh-Gordon model with even N, we have the following for- 
mulae for the operators A(u) and D(it) (recall that n = for the modular magnet and k = i/2 + bs 
for the sinh-Gordon model with even N) 

N-1 

A(u) = 2 N e" 7rKN sinhTrb(u + p N + p N ) *[[ smhnb{u - y fc ) (127) 

k=l 

\- t r sinh7r6(u - y;) , y,N 
fc=i ;#fc u /iy 

N-1 

D(w) = 2 N e~™ N sinh7r6(M + p N + p K ) sinb.7r6(u - y fe ) (128) 

fc=i 

+eii ^y^'t wx. 

^ smh7r6(y fe -y ; ) v ;/ 

where now p N = 2fe=i Yk't notice that [p N ,p N ] = and [p N ,y ] = -%t- Furthermore, compar- 
ision with ( IC. 13b and ( IC. 15b shows that p N = — J2k=i P fe ^ or tne m °dular magnet and pn = 
SfcLi( — ) k Pk for the sinh-Gordon model. 



6. Concluding remarks — outlook 
6.1 On the Baxter equations 

Summarizing the results of Sections @] and |3J we arrive at the main result of the present article. We 
will formulate it only for the case of the sinh-Gordon model with N odd for which our analysis is 
most complete, but from our previous discussions and remarks it seems clear that very similar results 
should hold in the other cases as well. 

Main result: A function t(u) is an eigenvalue of the transfer-matrix T SGo (u) if and only if there 
exists a function qt(u) which satisfies the following conditions 

(i) qt(u) is meromorphic in C, with poles of maximal order N in T_ s U T s , 

{exp ( + i7rN cr u — i|Nm 2 ) for |u| — > oo, |arg(u)|<§, 
exp ( — iw'Nau — i^Nw 2 ) for \u\ — ► oo, |arg(u)|>-|. 

(iii) t(u)q t {u) = (a(u)) q t (u - ib) + (d(u)) N q t (u + ib) , 
where d(u) = fi(-tt) = 1 + (l^) 2 e -*«>(2tH-») > 

(iv) qt{u) satisfies the following quantum Wronskian relation 

q t (u + iS + ) q t (u - iS + ) - q t {u + id-) q t {u - i8J) = Wjss(u) , 
where W N (u)= e-™ N ( u2+a2 )(D a {u)y N . 
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The corresponding eigenstate ^S>t in the SOV representation defined in Section \5\ is represented as 



We have therefore succeeded in reformulating the spectral problem for T(u) as the problem to de- 
termine the set 6 of solutions to the Baxter equation which possess the properties (i)-(iv) above. 

It should be observed that conditions (i)-(iii) already constrain the possible functions q t (u) rather 
strongly. Let us consider 

Q t (u) = (r 6 (f - i{u + s))r 6 (f - i{u - s)))~ g t (u) , (130) 

where Tb(x) = T2(x\b~ 1 , b), with T2 being the Barnes Double Gamma function defined in Ap- 
pendix lATl The function Qt(u) will then have the properties 



(a) Qt(u) is entire analytic of order 2 in C, 

(b) t{u)Q t {u) = (A{u)) N Q t (u-ib) + (D(u)) N Q t (u + ib) 



(131) 



The explicit form of the coefficients A(u) and D(u) can easily be figured out with the help of the 
formula d69|i and the functional relations dA.4l). 



Property (a) combined with the Hadamard factorization theorem (see e.g. O) imply that Qt(u) 
can be represented by a product representation of the form 



DO 



Qt(u) = e r ^H'[l--\ , (132) 

fe=i ^ Uk ' 

where the prime indicates the canonical Weierstrass regularization of the infinite product. The func- 
tion r(u) in the prefactor is a second order polynomial which can be worked out explicitly. The 
Baxter equation ( |57| | then implies that the zeros Uk must satisfy an infinite set of equations, 

x _ {A{u h )r Q t (u k -ib) 

which may be regarded as a generalization of the Bethe ansatz equations. However, as it stands it 
is not quite clear if these equations represent an efficient starting point for the investigation of the 
spectrum of our models. 

The quantum Wronskian relation (iv) encodes remarkable additional information which can not 
easily be extracted from the conditions (i)-(iii) above. We plan to discuss its implications elsewhere. 



6.2 Continuum limit 

It is certainly interesting to discuss the consistency of our results with existing results and conjectures 
on the sinh-Gordon model in continuous space-time. Let us therefore now show that our findings 
are consistent with Lukyanov's remarkable conjecture IILul on the ground state wave-function for 
the sinh-Gordon model in the SOV representation. 
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Recall from the Subsection 13.31 that we are interested in the limit N — ► oo, A — * 0, s — ► oo such 
that m = ^e~^ hs and R = NA/27T are kept finite in the limit. We are interested in the limiting 
behavior of the Baxter equation and of its solutions. Let us first note that the poles of q t (it) move 
out to infinity when s — ► oo. Also note that according to property ( 11291 -ii) rapid decay is found only 
within the strip 5 = {u€ C; |Im(tt)| < Q/2}. By noting that mA = <D(l/N) in the limit under 
consideration one sees that the coefficients a(u), d(u) in the Baxter equation ( 11291 -iii) become unity 
when N — ► oo. Most importantly, let us finally observe that the right hand side of the Wronskian 
relation ( 1129J -iv) approaches a constant for N — > oo. 

Our results therefore strongly suggest the following conjecture on the conditions which characterize 
the spectrum of the continuum sinh-Gordon model in the SOV representation: 

(i) qt(u) is entire analytic, 

(ii) qt(u) decays rapidly for |Re(u)| — > oo, ueS, 

(iii) qt (u) satisfies a difference equation of the form 
t(u)q t (u) = q t {u - ib) + q t (u + ib), (134) 

where t(u) is periodic under u — ► u + ib , 

(iv) q t (u) satisfies the following quantum Wronskian relation 
q t (u + i6 + )q t (u — i5 + ) — q t (u + i6-)q t (u — iS-) = 1. 

Our next aim will be to show that, by adding one supplementary condition, one gets a complete 
characterization of the function q (u) which was proposed in MLul to describe the ground state for 
the continuum sinh-Gordon model in the SOV representation: 

(v) qo(u) is nonvanishing within S. (135) 

We claim that the solution to conditions (i)-(v) is essentially unique and given by the formula 

mRcosh^u f dv log(l + Y(v)) 

(136) 



loggo(u) 2 s . n K b < J 2 Q cosh%(u-v)' 

which expresses qo(u) in terms of the solution Y(u) to the nonlinear integral equation 

/ j 
— s{u-v) iog(i + y(«)), 



(137) 



where the kernel S(u — v) is explicitly given as follows 

S(u) 



2 sin ^b cosh 



sinh ^(it + ib) sinh ^(u — ib) 



(138) 



These equations form the basis for the calculation of the ground state energy |Za| and other local 
conserved quantities of the continuum sinh-Gordon model [Lu| within the thermodynamic Bethe 
ansatz framework. 
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For the reader's convenience we will present the outline of an argument which establishes the 
equivalence between (i)-(v) and ( 1136b . (1137b . Let us define an auxilliary function Y(u) by the 
formula 

1 + Y(u) = q Q (u + iS+) q (u - iS+) . (139) 

Assuming that qo (u) satisfies the properties (i), (ii), and (v), one can take the logarithm of J 1 39b 
and then solve the resulting difference equation by Fourier transform, which leads to the representa- 
tion d 1 36b . Re-inserting this representation into the Wronskian relation (iv) shows that Y(u) must 
satisfy dLTTt . 

A proof of Lukyanov's conjecture [Lu| therefore amounts to showing that qo(u) must satisfy the 
properties (i)-(v). We find it very encouraging that our study of the lattice sinh-Gordon model gave 
us strong support for the necessity of properties (i)-(iv). 

6.3 Connection with lattice Liouville model 

Relations between the compact XXZ chain, lattice sine-Gordon model, and the (imaginary field) 
Liouville model was investigated in |FT| from the view point of the QISM. Let us show that similar 
connections exist between the modular XXZ magent, lattice sinh-Gordon model, and the (real field) 
Liouville model. Following BFTI . we introduce the L-matrix 

tf(u) = e-* 7 ^ 3 L SG {u + C) e^ K ° 3 , (140) 

where £ is related to the representation parameter s and the lattice spacing A as = Ae 7rbs , 
and the operators p and x in (l33l are related to the discretized field and its conjugate momentum as 
follows 

27r6p = 2ttK-/3$, 47r&x = 27r6C + /?(5n-$), = b\/8n . (141) 

Comparison with ( f47b and d49l shows that the new variables defined by ( 1141b are related to these of 
the sinh-Gordon model via a canonical transformation, 

n = n SG , /?$ = /?$ SG + 2tt&c, e -* K = %. (142) 



In the £ — * +oo limit L-matrix (1140b turns into 



ef n - (l + A 2 e-0*")ef n » -iAe" b 



V ; C^+oo v ' y _2iA sinh(7r6M+ f* n ) e"* 11 ™ J 

It is natural to expect that this L-matrix describes some massless limit of the sinh-Gordon model. 
The corresponding U-matrix obtained according to formula d52l ), 

/ &U(x) _i e «-|*(s) \ 

U (U)= I -2ismh(u+mx)) -|n(x) )' tl44) 



4 This argument is inspired by the considerations in [Za] and a suggestion of F. Smirnov (private communication, see also 
| Sm2|). However, the key point of our argument, namely the origin of the quantum Wronskian relation (iv) seems to be new. 
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reproduces the Liouville equations of motion via the zero curvature equation (see flFTI for details in 
the case of sine-Gordon model and imaginary Liouville field). This observation suggests that (11431 ) 
is a suitable L-matrix for describing the quantum lattice Liouville model in the QISM framework. 

Although the limiting procedure in (11431 ) has not been mathematically rigorously developed yet 
(in particular, there is a subtle question of interchangibility of ( — > +oo limit with the classical 
limit), the results of the present article provide further support for the proposed connection between 
the sinh-Gordon and Liouville models. First, observe that the twist by e5 7rb ^ CT3 and the shift of 
the spectral parameter in ( 11401 ) do not change the corresponding auxiliary R-matrix d34j l and the 
fundamental R-matrix given in Proposition Q] Therefore, the local lattice density of the classical 
Hamiltonian corresponding to d 140b can be obtained by substituting ( 1142b into d46| i. Taking then the 
limit £ — > +oo, we obtain the following lattice Hamiltonian density (up to an additive constant) 

H^ +1 ee Hm = ± log(± cosh § (H„ + II n+1 ) + \ cosh § ($„ - $„ +1 ) (145) 

+ a! e -|(*„+$„ +1 ) ^ + e !(n„+n„ +l) cosh _ <j, n+1 ))) t 

which in the continuum limit ( f50b yields the Liouville Hamiltonian: 

Es^;« + i^ const + P&dtf + KWf + ie-^). (146) 

The second observation that we can make to support the proposed relationship between the lattice 
sinh-Gordon and Liouville models is the following. The transfer-matrix corresponding to ( 1140b is 
given by T^(u) = T SG (u+C). Therefore, Q±(u) = Qjffu+C) satisfy (cf. (|79l l) the Baxter equation 

T<(u) ■ Qi(u) = (a c (u)) N Qi(u - ib) + (d c (u)f Qi(u + ib) , 
where a c (u) = 1 + A 2 e*W"-*>) , d ( (u) = 1 + A 2 e -^(2«+4C+i&) _ 

Hence for T L (it) = limf_>+oo T^(u) and Q±(u) = lim^ +oc Q±(u) (the limits are meaningful if 
T^(ii) and Q£(u) are expressed in terms of <f> and n), we obtain the Baxter equation 

T» • Ql(u) = (1 + A 2 e ^u-ib)^ Q l (u _ ih) + Q l (u + . 6 j ; (14g) 

which coincide^) with the Baxter equation derived for the lattice Liouville model by a different 
method in HK1IIFKV1 . 

We finish by noting that in the continuum limit, N — > oo, A = 0(l/N), the coefficient (l + 
^2 e Tr6(2u-j6)~j N b ecomes unity. This suggests that the Baxter equation for the eigenvalue q\(u) of 
the Q-operator for the continuum Liouville model is 



i L (u) q\{u) = q\(u - ib) + q\(u + ib) , t L (u + ib' 1 ) = t L (u) , 
5 modulo notations, in particular, A = 1 in [Kl FKV | 



(149) 
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which coincides with that for the continuum sinh-Gordon model dl34l ~iii). 

However, it seems to be crucial to observe that the asymptotic properties of the function q\ (u) will 
certainly differ from those found in the case of the sinh-Gordon model. Indeed, for any model, the 
asymptotic properties of q t (u) are related to these of t(u). Comparing the structure of the L-matrices 
d33l and (1143b . we see that the transfer-matrix T L (it) corresponding to the latter L-matrix has asym- 
metric asymptotics for Re (it) — > ±00. This seems to be related to the fact that the sinh-Gordon 
potential cosh/3$(a;) is spatially symmetric while the Liouville potential e~@®( x > is asymmetric. 
As a consequence, we expect that the set of solutions to the Baxter equation ( 1149b describing the 
spectrum of the continuum Liouville model will be quite different from the set of solutions to the 
same Baxter equation which characterizes the spectrum of the continuum sinh-Gordon model. 

These observations seem to offer a key to the understanding of the relation between massive and 
massless theories from the point of view of their integrable structure. 



A. Special functions 
A.l Double Gamma function 

All the special functions that we have to deal with can be obtained from the Barnes Double Gamma 
function T2(a;|w 1 , u) z ) BBal . which may be defined by 

logT 2 (x\oj 1 ,uj 2 .) = I ^2 (x + n 1 uj 1 + n 1 uj z y f \ . (A.l) 
\ ni,ri z =0 / {=o 

The infinite sum in ( I A. 1 1 is defined by analytic continuation from its domain of convergence 
(Re(i) > 2) to the point of interest (t — 0). One may alternatively use the integral representa- 
tion 

C , , x 1 f e- xt \oe(-t) dt 
logr 2 (,|^ 2 ) = -B 2A x\^) + ^J c {1 _ e _ uit) l_^ t) j (A.2) 

where the contour C goes from +00 to +00 encircling counterclockwise, C is the Euler's constant 
and 

^ = - — — — - 1^—^ (A - 3) 

The integral is well defined if Re(w 1 ) > 0, Re(w 2 ) > 0, and Re(x) > 0. It satisfies the basic 
functional relations 

F 2 (x +Cl) 1 \u 1 ,u z ) rz-u* " 2 T 2 (a; + w 2 ki,^2) rr- " 1 rK 

— ' ' ZlT — ; -, — r = V 27T — ; -. (A.4) 



r 2 (x\ui,LJ z ) T(x/uj z y T 2 {x\lj 1 ,U! z ) T(x/uj 1 ) 

(T 2 (x\uj 1 ,uj z ))~ 1 is an entire analytic function of order 2 w.r.t. its variable x with simple zeros at 
x = —mijj± — nijj z , where m and n are non-negative integers. 
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A.2 Function Wb (x) 

In what follows we will be dealing with 



r 2 (§ - ix\b~\b) 



w b (x) = £ . (A.5) 

T 2 (% + ix\b-\b) 

In the strip |Im(a;)| < §, function Wf,(x) has the following integral representation 



^ = eXP i ~ J ft sinh^sinhl } ' (A ' 6) 

K+iO ' 

where the integration contour goes around the pole t = in the upper half-plane. This function 
is closely related (cf. Eq. iD.lb ) to the remarkable special function introduced under the name of 
quantum dilogarithm in [FK2| and studied in the context of quantum groups and integrable models 

in llF2llRu1l^lPT^ICTIOlBTllT2llV2l . 

Analytic continuation of Wb(x) to the entire complex plane is a meromorphic function with the 
following properties 

self-duality Wb(x) — w b -i (x) , (A. 7) 

Wb(x + ib^ 1 ) +1 

functional equation =— — — = 2 cosh(7r6 x) , (A. 8) 

Wb{x - §fe ±:L ) 

reflection property Wb(x) Wb(— x) = 1 , (A.9) 



complex conjugation Wb(x) — Wb(— x) , (A. 10) 

zeros/poles (w fc (a;)) ±1 = ±x £ {i^+nb+mb^ 1 ; n, m e Z- } , (A.ll) 

i 

residue Res Wb(x) = — , (A. 12) 

x=-i% 27T 

[ e -^(^ 2 +^(6 2 +b- 2 )) for ijpi ^ ^ |arg(a;)| < § , 
asymptotics w b (x) ~ < _ (A. 13) 

[ e+T^ +12(6 + b ~ )) for |x| -> 00, |arg(x)| > § . 

Notice that |io;,(x)| = 1 if x £ K. Therefore, Wb(O) is unitary if is a self-adjoint operator. 

The function Wb(x) allows us to define a whole class of new special functions. In Appendix iDl we 
will use in particular the following b-analogues of the hypergeometric functions defined by 

* r (U 1 ...U r ;Vi...V r ;x) = \ [ dr e™> f[ ^ + T \ , (A.14) 

% J f-J: bb{V k + r) 

where the special function Sb(x) is defined by 

Sb(x) = w b (ix - iQ) (A.15) 
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and has the properties 

self-duality S b (x) = S b -i(x) , (A.16) 

functional equation Sbix+b^) = 2 sinijrb^x) Sb(x) , (A. 17) 

reflection property Sb(x) Sb(Q — x) = 1 . (A. 18) 



A.3 Function D a (x) 

Let us also introduce another useful function 



D a {x) = ?f±°>. (A.19) 
Wb{x — a) 



Combining dA.6b with dA.221 >. we derive the integral representation 

„ , , { [ dt cos(2ia;) sin(2ai) 1 , A „ ON 

^=4i 2i sLismhf T (A ' 20) 

R+iO ' 

D a (x) is a meromorphic function with zeros at ±x e T_ Q and poles at ±x € T Q , where the set 
T Q is defined in (IB.lt . The function D a (x) is self-dual in 6 (but we will omit this index) and has 
the following properties 

e - , D a (x + ^b ±1 ) coshnb^ix + a) 

functional equation — — ITZTT = — r — txtt r, (A.21) 

D Q (a; - |b ±:L ) coshnb ±L (x ~ a) 

x-parity D a (x) — D a {—x) , (A. 22) 

reflection property D a {x) D_ a (x) = 1 , (A.23) 



complex conjugation D a (x) = D- & (x) , (A. 24) 



g-Ano* for | x | ^ 00) (arg^)! < | 



x-asymptotics L> Q (x) - <j +2maxc ' !"T'J (A.25) 

[ e+^" Q2; for | act) -> oo, |arg(a;)| > | , 

f e- iw (- 2 +- 2 + A(6 2 +b- 2 )) if | a) _ 0O) | arg(a) | < n ( 

a-asymptotics D Q (x) ~ <^ _ (A.26) 

I e +Mr(a: +« +T2(& +6 )) if | a | _> oo ; |arg(a)| > § . 

Also, the following identity is obvious from the definition dA. 19b 

Djx) Dn(v) = D a+ , r -« ( x +»+ a -P ) D a+P -^ v ( x+ K a+p ) ■ (A.27) 

Notice that |D a (a;)| = 1 if a € M and lelonejl. Therefore, D a (0) is unitary if a e K and 
is a self-adjoint or anti-self-adjoint operator. 
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A.4 Integral identities for D a {x) 

Here we will give some integral identities involving products of D-functions. These identities can 
be regarded as summation formulae for the b-hypergeometric functions $ r introduced in (IA. 14t . 

Let us denote a* = — | Q — a and introduce the function 

A{a\,ot2, ■■■) = Wb(ai-al)wb(a2-at2) ■■■ ■ (A.28) 
Notice that (a*)* = a and hence A(a*, 0%, , . . .) A{a,\, 0,2, ,...)= 1. 

Lemma 15 in MPT2II and Eqs. (26)-(27) in BFKV1 can be rewritten as the following property of the 
function D a (x) under the Fourier transform: 

f dx e 2mxy D a (x)= A{a) D a * (y) . (A.29) 

Taking into account that lim a ^o D a (x) = 1, we obtain from ( |A.291 > as a special case 

lim A(a*) D a (x) = S(x) . (A.30) 

Here S(x) on the r.h.s. is the Dirac delta-function and this relation should be understood in the sense 
of distributions. Indeed lim Q ^_^Q A(a*) = Wb(^Q) = and the l.h.s. of JA.301 > vanishes almost 
everywhere. On the other hand, the only double pole of D_±q(x) is at x = 0. 

Using ( |A.29t , it is easy to derive the following relation 

I dx e 2mzx D a (x - u) D p (x - v) = 

(A. 31) 

= A(a,P)e^ u +^ [ dye 2 ^ u -^D a .(y + z)D >(y-±). 

Choosing z = a* + (3*, we can use (|A.27t in order to rewrite the product of D's in the integrand on 
the r.h.s. as a single function, D a ± + p> (y + a ~@ ), and then apply ( |A.29t . This yields 

/ dxe 2vi! - a * +f3 *) x D a {x-u) Db{x-v) = 

Jr (A.32) 

= A(a, /3, a*+p) e ^(««*+^*) D a+0+ j Q (u - v) . 

In the case a* = —(3*, Eq. (1A.301 I can be used and we conclude that 

/ dx D a (x - u) D p (x -v) = A(a, /3) 5(u - v) (A.33) 

JR 

holds in the sense of distributions provided that a + (3 = —iQ. 
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Lemma 3. The identities 

dxD a {x—u)Dp{x—v)Dy(x—w) (A. 34) 

= A(a, j3, 7) D a * (w—v) Dp* (u—w) D 7 * (v—u) , 
dx D a {x—u) Dp(x—v) D 1 {x—w) D UJ (x—z) (A. 35) 

D a+ p + i_ Q (u-v) r 

A(a, /3, 7, w) — - — - f - / dx D a * (x— v) Dp* (x— u) D 1 * (x— z) D u * (x— w) 



D a+p+iQ {w-z) 

are valid provided that a + (3 + 7 = — iQ in HA.34}) , and a + P + j + u> = —iQ in ( IA.35I ). 

Proof. Relation (IA.341 follows straightforwardly from Eq. (19) in HK2I . where function ipb(x) is our 
g b (e 27rbx ) (cf. Eq. dDJl i). Also, in other notations, relation ( lA34b is Eq. (1 1) in lTV2l . 

Eq. (IA.34b provides two expressions for a function which we denote as I(u, v, w; a, (3, 7). In order 
to prove ( IA.351 we multiply two copies of ( IA.341 and compute the following integral 

dt I(t, u, v; a, (3) I(t, w, z; /i, 7, uj) , (A.36) 

where v + a + [3 = fi + j + uj = —iQ and we impose an additional condition v + p, = —iQ. Then 
the l.h.s. of ( IA.35t is obtained if we substitute for the /'s the expressions on the l.h.s. of ( I A. 341 ) and 
use relation ( IA.33b . The r.h.s. of ( IA.35t is obtained directly from ( IA.36t if we substitute for the /'s 
the expressions on the r.h.s. of ( 1A.341 > and use also that ji* = —v*. □ 



B. Positivity versus self-duality of the representations V s 

The representations V s are distinguished by the property that the operators ir s (u), u G {E, F, K} are 
positive self-adjoint. We are now going to show that this property is closely related to the remarkable 
self-duality of these representations under b — > b^ 1 which has such profound consequences for the 
physics of our models. 

To begin with, let us remark that there exists a linear basis B q (sl(2, M.)) for U q (sl(2, R)) such that 
all elements u of B q (st(2,R)) are realized by positive operators n s (u). Such a basis is, e.g., given 
by the monomials 

q +n ^C l E m K n represented by (K s )t (E s ) m (K,)*, 

mn , n „ /, TO, 71 € Z, Z, TO > 0. 

q-—C l F m K n represented by (K s ) * (F s ) m (K s ) = , 

The elements of U q (sl(2, R)) are clearly realized by unbounded operators on i 2 (R). It is therefore 
useful to consider suitable subspaces T s C L 2 (R) of test-functions on which all operators 7r s (it), 
u £ U q (sl(2, R)) are well-defined. In order to describe a canonical choice for T s let us represent the 
elements of T s by functions f(k) such that p acts as (p/)(fc) = kf(k). 
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Definition 5. Let T s be the space of functions f(k) which satisfy e a ' fc ' / £ L 2 (R) for all a > 0, ana" 
which have an analytic continuation to C \ (T s U T_ s ), where 



T s = |s + i(f + n& + mb~ 1 ), n,meZ^ | 
T s = |s-i(f +nb + m^ 1 ), n,raeZ^ } 

On the spaces T, the action of U q (sl(2, R)) is given by 



(B.l) 



E. /(*)=[§ +w-ifc] 6 /(A : + i6), ^ f; _ M 
F s f(k)=[% + is + ik] b f(k-ib), 



K s f(k) = e-* bk f(k), (B.2) 



where \x] b = . 

The distinguished role of the space T s is explained by the following result, which shows that the 
space T s is canonically associated to the representation w s : 

Lemma 4. T s is the largest space on which all tt s (u), u € U q (sl(2, R)) are well-defined, i.e., 

%= p| V u , (B.3) 

uG8,(6l(2,R)) 

where T> u is the domain of the unbounded operator ir s (u), u G B q (sl(2, R)). The spaces T s , s G R 
are Frechet spaces with topology defined by the family of seminorms 

|| / || w = sup | (7r.(tx)/)(fc) |, «6B,(«[(2,R)). (B.4) 



Proof. It is easy to check that ir s {u)f S T, for all f £ T s . In order to show that the conditions in 
the definition of T s are all necessary, let us first observe that e a ' fe '/ € L 2 (M) is clearly necessary 
for K™/ to be well-defined for all n £ Z. In order to determine the conditions on / for E™/ to be 
well-defined, let us consider the unitary operator U s = Wb(p — s), where the special function Wb(x) 
and its properties are described in Appendix lA.2l We then have 

U s • E™ • U; 1 = e 27rnhx . 

The intersection of the domains of e 2mbx for all n £ N consists of functions g(k) that are analytic in 
the upper half plane H + , see, e.g., (Sj Lemma 1]. The corresponding functions f(k) = (Uj 1 g)(k) = 
(wb(k — s))~ 1 g(k) may have poles in T s . Similar arguments applied to F" allow us to complete the 
proof of the first statement in LemmaH 

In order to verify the second statement, we mainly have to show that the space T s is complete w.r.t. 
the topology defined by the seminorms ( IB. 4b . This follows from ( IB. 31 > together with the observation 
that the self-adjoint operators 7r s (u), u £ B 9 (fi[(2, R)) are closed on T> u . □ 
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We regard Lemma|4]as the key to the mathematical understanding of the duality b — > b^ 1 of our 
representations V s . Indeed, let us introduce the operators E s , F s , K s , obtained by replacing b — ► b^ 1 
in ©. These operators generate a representation V s of Uq(sl(2,M.)), q = e i7rb on the same 
space T s . The space T s is associated to the representation V s as canonically as it is associated to V 8 - 
Moreover, it is easy to see that the representation V s commute^ with V s on T s . It is therefore natural 
to regard T s as the natural space on which a representation of the modular double U q (sl(2, R)) (g) 
Uq(sl{2,R)) E3] is realized. 

Another way to make the self-duality of the representations Vs transparent uses the rescaled gener- 
ators introduced in (fTTT i. These generators and their counterparts e s , f s , k s , obtained by replacing b 
with b^ 1 are related as [BT| 

(e s )i=(e s )\ (f s )i = (f s ) b , (k s )i = (k s ) 6 . (B.5) 

These observations express quite clearly that the representations of the two halves of the modular 
double, U q (sl(2, R)) and Uq(sl(2, R)), are related to each other like the two sides of the same coin. 

C. Structure of the monodromy matrix 

This appendix is devoted to the derivation of some simple, but important structural properties of the 
monodromy matrix M(it), 

M(u) = ( An( ; U \ Bn( ;\) EE L N (u) L 2 (u) ■ Ll {u) . (C.l) 

C.l Expansions in the spectral parameter 

Introduce the following notations for n £ N, 

,_ In — 1 if n - odd; In if n - odd; 

M = { > W = { • (C2) 

In if n - even; I n — 1 if n - even; 

Lemma 5. The elements An (it), Bn(w), and Dn(w) ofM(u) have the following form: 

x 



A^ xz (u) = e^ bu {-) m e- 2mlTbu A*^ , (C.3) 

771 = 

N-1 

Bxxz/„,\ ■ „N7rfm \ * / \m „ — Imixbu dxxz tr^ a\ 

n W — ie / jK-) e B N,m' ( c - 4) 

771 = 

N 

D N xz (u) = e Nnbu (-) m e~ 2m * hu D N X ^ , (C.5) 



771 = 



6 Commutativity of V s and Vs only holds on the dense domain T s but not in the usual sense of commutativity of spectral 
projections ! 
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TNI 

a^ g (u) = *r*n e T»(fNi«-N«) J2 (-)™ e - 2m7Ibu a s n % , (c.6) 

m=0 
LNj 

B^ G (u) = — iLNJ e 7rb(LNj„-N s ) ^ ^m^m^&u gSG m > (CJ) 

m=0 
TNI 

d^ g (u) = *r*n e -KrNi«-N S) ^ ( _ )me -2m^« D so m ) (C 8) 

m=0 

where An,™, Bn, to . an d Dn,™ are positive self-adjoint operators. 

Proof. Let us consider the case of the XXZ chain, the other case being very similar. The definition 
of the monodromy matrix M N (u) yields the following recursion relations 



A N X » 


= (e w6u k N - 




h «e t N l- N _i(U,), 


(C9) 


B N X » 


= (e^kN- 






(CIO) 


Q x » 






-ie-'^eNA^Cu), 


(Cll) 




- (e^k- 1 


-e-^k N )D™( U )- 


f » e -^e N B™( U ), 


(C12) 



where k]y = k a <g> 1 (g> . . . 1 etc. Using these recursion relations one may inductively show that the 
operators An.,„, Bn,™ and Dn, to are linear combinations of monomials of the form 

uf) ® . . . ® u5e{e a ,f a ,k a ,k7 1 } 

with positive integer coefficients. It remains to note that an operator which is the sum of positive self- 
adjoint operators will be self-adjoint on the intersection of the domains of the individual summands. 
These obserations reduce our claim to the self-adjointness and positivity of e a , f a , k a . □ 

Lemma 6. The leading terms o/An(m), Bn(u), and Dn(u) at e* hu — > ±oo are given by 

B^ z = A^- 1 ) f , A N * Z = 0$$ = A**" 1 ) k , 

B^_! = A^" 1 ) f , A^ = D N * Z = A^" 1 ) k- 1 , 



(C.13) 



where A^ is the n-fold co-product defined via A( n+1 ) = (A^™) ® id) o A, with A^ = id and 
A^ 1 ) = A, and A^ is defined analogously for the opposite co-product A^ = A'. 

[ b^ g = 9 odd (a( n -d k- 1 ) , b- n = e odd (a^- 1 ) k) , 

N-odd: I A^ = o dd (A^ N_1 ) f) , A^Vi = fodd(A( N - 1 )f), (C.14) 
{ D^ = odd (A^- 1 ) e) , D^ N _! = 6 odd (A^- 1 ) e) , 

f B^ G = 9 cvcn (A( N -D f) , B^ J G N _ 1 = 9 cvcn (A( N -D f) , 
n-even: < A|j G = 9 cvcn (A( N_1 ) k) , Af, G N = 9 cvcn (A^ k" 1 ) , (C.15) 
I - 9 cvcn (A( N -D k- 1 ) , D^ G N = 9 cvca (A( N -D k) , 
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where 6> dd = #n-i o . . .630 6\ and t 
( 1571 ) at odd/even sites. 



0N • • • @4 O2 are compositions of the automorphism 



Proof. Eqs. ( IC. 1 3b follow easily from the decomposition ( fT4b of the corresponding L-matrix. 
Eqs. ( IC.15I ) are obtained by analogous consideration if L SG (u) is replaced with o\ L'(u) (see (l4TT i) 
and formula d42l > is used. In order to apply this approach in the N odd case, one has to multiply 
the monodromy matrix with an extra o\ from the right (which leads to the interchange A <-> B and 
C <-> D). □ 



C.2 Quantum determinant 

Let us discuss connection between the so called quantum determinant and coefficients a(u) and d(u) 
which arise in the Baxter equation d62l i. Since for the modular magnet we use R-matrix ( fTBI l which 
is not symmetric, the corresponding quantum determinant will differ from the "standard" formula 
applicable, e.g., for the sinh-Gordon model. Therefore we commence by deriving the required 
expression. 

Lemma 7. Let L{u) be an L-matrix, satisfying the relation ( 1721 ) with the auxiliary R-matrix of the 
form 



( sinh7r6(u + ib) 



R(u;0 = 



sinh irbu i sin irb 2 e"^" 
i sin7r& 2 e _7rb ^" sinh7r&u 



V 



sinh7r6(u + ib) J 



(C.16) 



and let M (u) be the corresponding monodromy matrix defined by iC.lj . The following element 
(quantum determinant) 



det g M(u) = A(u) D(u - ib) - B(u) C(u - ib) 
is central, i.e., [M(v),det q M(u)] = 0, andean be written as follows 

det q M(u) = (det g L(u)) N , 



(C17) 



(C.18) 



where the quantum determinant of L(u) is defined by the same formula \C.17\ (with M(u) replaced 
by L(u)). 

Proof. Existence of the quantum determinant is due to the degeneration of the auxiliary R-matrix 



I 



R(—ib;^) = isinirb 2 



-1 



\ 



(C19) 
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It is interesting to notice that this matrix is proportional (in the standard basis) to the one-dimensional 
projector onto the spin representation in the tensor square of spin ~ representations of 
U q t(su(2)). 

For R(u; 0) the statement of the Lemma is well-known (see, e.g., II KB ID . In the generic case £ ^ 0, 
one can observe that the gauge transformation 

L(u)=g- 1 L(u)g u , M(u) - g' 1 M(u) g u , g u = e^ h ^ (C.20) 

yields L-matrix and monodromy matrix which satisfy the relation ( fT2b with the auxiliary R-matrix 
R(u; 0). Therefore 

det, M(u) = A(it) D(u - ib) - B(u) C(u - ib) (C.21) 

commutes with entries of M(u) and hence with entries of M(it). Using ( 1C.20I ) in order to rewrite 
( 1C.2U in terms of entries of M (u), we obtain ( IC. 17b . 

In order to prove (IC 18b it suffices to observe that 

det 9 M(u) ■ R 12 (-ib;£) = M^u) M 2 (u - ib) Ri 2 (-ib; £) , (C.22) 

which yields also three different expressions equivalent to dC. 17b if we take into account the relation 
Ri2{-ib;C) Mi(u - ib) M 2 (m) = Mi(u) M 2 (u - ib) -R12 which is a particular case of (fT2l . 

Now if M'(m), M"(u) satisfy (fT2l) with the same R-matrix and their entries commute, then we have 

det, (M'(u) M"(u)) ■ R 12 (-ib; £) = M[(u) M'{(u) M' 2 {u - ib) M 2 '(u - ib) R 12 (-ib; 

= M[{u) M' 2 (u - ib) M'{(u) M' 2 \u - ib) Ri 2 (~ib; = dct q M'(u) ■ dct q M"(u) . (C.23) 

Whence JC. 1 8b follows immediately. □ 

For the models that we consider, Lemma|7]yields 

det 9 M xxz (w) = A(u) D(u - ib) - q' 1 B(u) C(u - ib) (C.24) 

= (4 cosh7r6(s + 11- i|) cosh7r6(s -u + «|)) N , (C.25) 
det 9 M SG (w) = A(u) D(u - ib) - B(u) C(u - ib) (C.26) 

= (4e^ bs cosh7r6(s + u - i|) cosh7r6(s - u+ i|)) N . (C.27) 

Proving the Proposition [5] in Appendix E, we will deal with monodromy matrices M(u) obtained 
from M(u) by a gauge transformation with a matrix G(u). Notice that such M(u) satisfies the 
exchange relation ( flZb with an R-matrix 

R 12 (u, «; = Gi (u) G 2 (v) R(u-v; G^(u) G 2 \v) (C.28) 

which is not of the form ( IC. 16b . Therefore, for general G(u), Lemma [7] does not apply to M(u). 
Nevertheless, there exists a class of gauge transformations which preserve the quantum determinant 
in the following sense. 
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Lemma 8. Let M(u) satisfy ( 1721 ) with the auxiliary R— matrix R(u\ £) of the form AC.16h Let M(u) = 
( <5(u) ) ^ e !fs S au 8 e transform defined by 

M(u) = G(u)-M(u)-G- 1 (u), G(u) = ( 1 p ° 1 ), p = e^ 11 p , (C.29) 
where po is a c-number, which does not depend on u. Then 

det 9 M(u) = k(u) D(u - ib) - B(u) C(u - i&) = det g M(u) , (C.30) 
where the r.h.s. is defined according to Lemma^\ 

Since entries of M (u) are linear combinations of entries of M (it), it follows that det g M (u) is central, 
i.e., [M(u), det g M(u)] = 0. Thus, the l.h.s. defines the quantum determinant corresponding to the 
R-matrix dC28l 

Proof. Although R(u — v: £) and R(u, v; £) are in general not equal, they coincide at u — v = —ib. 
Indeed, using the explicit form of G(u), it is easy to check that 

Gi(u)G 2 (u-ib)R(-ib;£) = R(-ib;£). (C.31) 

Observing also that ( IC.22I ) holds for det g M(u) as well, we derive 

det q M{u)-R 12 (-ib;Z) = Mi(u) M 2 {u - ib) R 12 (-ib;£) 

= Gi(u) Mi(u) G^{u) G 2 {u - ib) M 2 (u - ib) G^{u - ib) R 12 {-ib; £) 

= Gi(u) G 2 (u - ib) Mi(«) M 2 (u - ib) 1 (« - ib) G^(u) R 12 (-ib;£) 

= Gi{u) G 2 {u - ib) Mi(u) M 2 (u - ib) Ri 2 (-ib; 

= Gi(u) G 2 (u-ib) Ri 2 (-ib;C) -det g M(u) = Ri 2 (-ib;Q ■ det q M(u) , 

which proves the assertion of the Lemma. □ 

It is important that the gauge transformations used in the proof of the Proposition[5]belong to the class 
of gauge transformations described in Lemma[8j they correspond to £ = 1 and £ = 0, respectively. 
This fact allows us to relate the quantum determinants of the models in question and the coefficients 
of the corresponding Baxter equations. A quick inspection of the proof of Proposition[5]shows that 

(a(u) d(u - ib)) N = det g M (u) , (C.32) 

where the r.h.s. is given by ( IC.25I ) and ( IC.27I ). respectively. 

D. Construction of the fundamental R-operator R(u) 



For the proof of Theorem[T]we will need the following material from 1PT2I IBTll . 
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D.l Clebsch-Gordan maps 

Let the space M. be defined by the direct integral 

M = / dsV s . (D.l) 



+ 



Realizing elements of T s as functions f(k) leads us to represent the elements of M. by families 
of functions / = ( f s ; s £ R + ), where f s = f s (k) £ T s for all s £ R + . We shall define the 
multiplication operator s by 

sf = (sf 3 ; s£R+). (D.2) 

To any family ( S ; s £ R + ) of operators on T s we may then associate an operator S on A4 in the 
obvious manner. We have the corresponding canonical action of U q (si(2, R)) on M. via 

ir s {X)f = (n s (X)f s ; sGR+), V X G W g (s[(2,R)). (D.3) 

The Clebsch-Gordan maps C SlSl were defined in IPTll lPT2l as a family of operators 

C SlSl :? Sl ®P Sl ^M. (D.4) 

The Clebsch-Gordan maps C SzSl intertwine the action of U q (sl(2, R)) on 7- > Sz (giT 3 ^ with the canon- 
ical action on A4 in the sense that 

C 5iSl • (tt 52 ® tt s JA(X) = 7T S (X) • C SlSl VXeW,(5l(2,K)). (D.5) 

D.2 R-operator and braiding 

Let us introduce 

R s , Sl = g H *-® H «i 5h (e Si ® f.J g H -»® H -x . (D.6) 

Here the anti-self-adjoint operator H s is defined by k s = g Hs , and gb(x) is the non-compact quan- 
tum dilogarithm related to the special function Wb(x), defined in AppendixlAl via 

g b (e 27rbx ) = e % ib2+b ~ 2)+ ^ x2 Wb(-x) . (D.7) 

The R-operator R SiSl satisfies the following relations [BT| 

R SiSl • (tt S2 (g> tt Si )A(X) = (tt S2 ® 7r Si )A'(X) • R S2Si , (D.8) 
R SiSl (e Sl ® k- 1 ) = (e Si ® k,J R SlSl , (D.9) 

R s , Sl (k Si ® f.J = (k- 1 ® f 8l )RM, , (D.10) 

where X £ U q (sl(2, R)) and A' stands for the opposite co-product. 

The braiding operator B : V Sz ® V Sl — > V Sl ® Ps 2 is defined by B SiSl = P R SiSl > where P is the 
operator that permutes the two tensor factors. In what follows we will need the following statement. 
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Proposition 3 (Theorem 6 in [BT|). The braiding operator is diagonalized by the Clebsch-Gordan 
maps in the following sense: 

C SlSi • B SiSi = 0/ iSi • C S2Sl , (D.ll) 
where f2 s s w operator on M associated via ( IP. 2D to the scalar function 

■l 2 i 2 2 i Q 2 \ 

In the particular case, s z = s 1; one can regard the permutation P as an endomorphism of V Sl <8> V s ± 
and then Proposition^ allows to relate it to the R-operator: 

P = R SlSl e -™(^+2^ 2 ) . ( D .i3) 

D.3 Proof of Theorem Q] 

We will consider the more general R-operator defined onV Sz <8 V Sl by the formula 

R 8i sAu) = R SiSl e +"( 5 "- s2 1 - s '-T) Du ( Sai ) (D.14) 

where the special function D a (x) is defined by Eq. flA.19t in Appendix [A] and s 21 is the unique 
positive self-adjoint operator such that 

4cosh 2 ?r&s 21 = (tt Sz <g)7r s JA(C) . (D.15) 

It follows easily from relation ( ID. 13b that R(u) = R ss (u) coincides with the fundamental R-operator 
defined earlier in (|20T >: 

R{u) = P w b (u + s) w b (u - s) = P D u (s) . (D.16) 

To begin the proof of TheoremQ] which is somewhat more involved than the proofs of the analogous 
results in the case of highest weight representations fjj] IF1I IBD1 . let us observe that Eq. $1% is 
equivalent to the following set of equations 

R SiSl (u) A(L ± ) = A'(£ ± ) R S2Sl (u) , (D.17) 
R SiSl (w)^(u)=r(£ ij (-u))R S2Sl («) ! i,i = l,2, (D.18) 

where were introduced in ([T3l-(fl4li. r is the flip operation: T(a S2 ®6 Sl ) = b Sz ®a Sl , and the 
operators £y (it) are given by 

* u (u) = e* bu e Si <g> f Sl + e* bu k" 1 ® k Sl + k Si ® k" 1 , (D.19) 

M«) = e* bu k Si ® f Sl + e^ 6 " f Si ® k" 1 , (D.20) 

M«) = e" bu e Sl ® k- 1 + e^ bu k Si ® e Sl , (D.21) 

£ zz (u) = e^ 6 " f Si ® e Sl + e^ 6 " k- 1 ® k Sl + e Tta k Si ® k" 1 . (D.22) 
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The verification of the relations ( ID. 171 ) is easy. Introducing 

r s >M = Kk-KsA u )> (D.23) 
it follows from ( ID. 8b that ( ID. 17b is equivalent to the system of equations 

r S2Sl («) A(X) = A(X) r SiSl (u) . (D.24) 
Validity of the relation dD~24i i follows from [s, A(X)] = 0. 

The verification of the relations (ID. 18b is somewhat harder. To begin with, let us observe that it 
suffices to verify relation (ID. 20b . say. 

Lemma 9. Validity of \D.18$ for (i, j) = (1,2) implies that R s , Sl (u) satisfies the three remaining 
relations in AD.18\l as well. 

Proof. Equivalence of the (1, 2) and (2, 1) equations in ( ID. 18b can be established by invoking the 
automorphism 6. defined in (l37b. Introduce t 9 = (6®6)ot. It is easy to see that r e (R SiSl ) = Rs^ 
and r 6) (A(C)) = A(C). Therefore, T 9 (R SiSl (it)) = R SzSl (u). The claimed equivalence of the 
(1, 2) and (2, 1) equations in (ID. 18b follows now by observing that r e (£ 12 (u)) = £ zl (u). 

In order to prove that relation ( ID. 18b for (i, j) = (1,1) follows from the validity of relations ( ID. 18b 
for (i, j) — (1, 2), (j, j) — (2, 1) and ( ID. 17b , let us consider the following object: 

X(u) = q(k 2 Sz ® k 2 s± - k- 2 ® k- 2 ) £ lx (u) + k Sz ® k Si (e-* bu C Sl <g> 1 + e" bu l ® C.J 

+ k" 1 ® k-^e^C^ ® 1 + e-^l ® C a J . (D.25) 

The operator (k 2 ^ ® k 2 ^ — k;r 2 ® k~ 2 ) does not have a normalizable zero mode. Validity of relation 
( ID. 18b for = (1,1) therefore follows from 

R SzSl (u) X(u) = r(X(-«)) R SiSl («) . (D.26) 

Validity of this relation follows from the observation that X(u) can be represented in the form 

X(«) = (k S2 <g> k Si )e iz (u) (7r Si ®7r s JA(e) - (k" 1 ® k" 1 )^^) (7r Sz ®7r s JA(/) . (D.27) 

Exchanging k «-* k _1 in ( ID. 27b , we can derive relation (ID. 18b for = (2, 2) in a completely 
analogous way. □ 

It remains to verify relation ( ID. 18b for = (1,2). This relation may be rewritten in terms of the 
operator r SzS^ (it) defined in dD.23b by using the relation ( ID. 9b . We conclude that validity of ( ID. 181 ) 
follows from the validity of 

r S2Sl (t») (e^Oi + O a ) = (e 2 ^ u 0, + 3 ) r SiSi («) , (D.28) 
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where we have introduced the convenient abbreviations 

Oi = e s , (g) k" 1 , 

n , " 3 ^R- 1 Si (k s ->e s JR SiSi . (D.29) 

In order to prove that r SiSl (u) satisfies ( ID. 281 ) it will be convenient to use the Clebsch-Gordan 
maps C SlSl introduced in Subsection lD.il Let us observe that dD.241 > implies that the r-matrix is 
diagonalized by the Clebsch-Gordan maps: 

C SiSi • r Sl5 >) = < s >) • C Sl5l , (D.30) 



where r| a (u) is the operator on M. associated to the scalar function r/ a (u) via ( ID. 21 ). 

In order to further evaluate equation dD.28b we will need to describe the images Og, £=1, 2, 3 of the 
operators Oi under the maps C SlSl which are defined by C SiSl • 0^ = 0^ • C S2 Sl . 

Proposition 4. The operators Og can be represented as follows: If g = (g s ;s £ K + ) G M, then 
6 e g ee ( {6 e g) s ; s € R+ ), where 

l 

(6 e g) s (k) ee A Z lSl ( fc ) iy s9s(k+ib) , (D.31) 

i/=-i 

where T"g s — g s +iub- The coefficients A^^ lSl (k) are symmetric under exchange of s z and s±, 

A^^(k)=A^ s >(k) (D.32) 
and are otherwise related to each other by 

2 ;s 1;s (D.33) 

A^' SiSl (k) ~ e - 27Tbl/s - i7Tb2 v 2 A^'< s t s i (fc) 



The proof of Proposition|H which is somewhat technical, is given in Subsection ( ID. 51 ). 

Proposition [4] together with equation ( ID.30I ) allows us to rewrite the defining relation ( ID. 281 ) as a 
commutation relation satisfied by the corresponding operators on A4. Applying ( ID. 281 ) to a func- 
tion g s and matching the coefficients in front of g s +i U {k+ib) in the resulting equation, we derive 
functional equations on r/ iSi (u), 

r£*» (e 2 ^A^(k)+A^(k)) 

= (e 2 ^" ^; ; s SlSl (*) + ^'V r ' 4x (u) . 

The case f = holds trivially. Taking into account ( ID. 321 ) and ( ID. 331 ). Eqs. ( ID.34I I for v — ±1 are 
equivalent to a single functional equation, 

(e 2 " bu + C) = (e 2 ^ + C 1 ) < Sl («) , (D.35) 
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where £ = e i7rb + 27r6s . In terms of f/ s (it) = e _l,rs (u) one may rewrite this functional 
relation as follows 

r^a, 2 (u) cosh7r6(s +u) = f 8 Sz tf (it) cosh7r&(s — u) . (D.36) 

Recalling that the special function Wb(x) satisfies (IA.8I 1. we conclude that equation (ID. 35b is solved 
by the expression for r/ which follows from eqn. ( ID. 141 ) via ( ID.23I I and ( ID. 30b . namely 

r U (u) = e >«*-*-*-& = e^ 2 -^-^) D u (s) , (D.37) 

- Wb(S — u) 

where D a (x) is defined in (|A.19t . 

Property JA.24I > implies that |? , / iSl (u)| = 1 if u, s E M. Since s 21 is self-adjoint, we infer that 
r s s jSi (u) is a unitary operator for «el On the other hand property (|A. 10b implies that gb{x) given 
by dD.7| i satisfies |5b(x)| = 1 for x G M + . Since e 2 ® in (|D.6t is positive self-adjoint, we infer 
that R SlSi is unitary. The unitarity of R SzSl (it) for u€l follows. □ 

Let us finally remark that relation jA.261 ) leads to the following asymptotics of R SlSl (it) 
R SlSl (u) ~ e -*'(»"+-) R SiSi for Re(n) -» +00 

(D.38) 

R S2Sl («) - e +^ u +-) R SiSl for Re(u) -> -00 , 
where R_ , = P ■ R7 1 , ■ P. 



D.4 Clebsch-Gordan and Racah-Wigner coefficients for V s 

In order to prove Proposition [4] we will need to describe the Clebsch-Gordan maps C SlSl more 
explicitly. For the following it will be convenient to use the variables a r — § + is r in order 

parameterize the representations V Sr , r = 0, 1, 2 The Clebsch-Gordan maps C SzSl can then be 

represented explicitly as an integral transformation 1PT2I BT | : 

/.(*) = J dk 2 dkr [Z I Z t ] /(**, fc), (D.39) 

The kernel which appears on the right hand side, the so-called b-Clebsch-Gordan kernel, was cal- 
culated in IPT2l lBTl . It is of the general form 

[tl I % Z } = S(k 3 -k 2 - h) [as I Z Z } (D.40) 
where the function [ «3 | Z Z ] i s given as 

Tri (1) (0) /, , \ 

Hfa *i] = g , ^ — ^$ 3 (Ri,R2,R 3 ;S u S 2 ,Sr,-a[°> 3 ). (D.41) 
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In equation ( ID. 401 ) we have used the special function <j> 3 (i?i, R 2 , R3; Si, S 2 , S3; x) defined in ( IA. 141 ) 
whose arguments have been chosen as follows: 

Ri = cti + iki Si = o^23 + ai + ik\ 
R 2 =a 2 - ik 2 s 2 = af 2 } 3 + a 2 - ik 2 
R3 = a°23 S3 = Q. 

Here and below we are using the following notations: 



(°) 11 r> 

a ijk — a t + a J + a k - Q- 



The precise relation to the b-Clebsch-Gordan coefficients from [BT| is as follows: 

r a 3 I a 2 cc\ I _ 013 r«3 I «2 «ll 

L fc 3 I k 2 fci J here — a 2 ai [ fc 3 | fc 2 fci J BT 



41 


= a j ~ 


V a k 


- Oti, 


(2) 
a ijk 


= OL k - 


f a% 


-Otj, 


(3) 

a ■ 

13k 


= OLi H 


-atj 


- a k , 



(D.42) 



with phases v^lai chosen as 

{Kiai) 1 = S b (a 3 + ai - a 2 )S b {a 3 + a 2 - oti). (D.43) 
The b-Racah-Wigner coefficients are then defined by the relation 

r a s I ai ao 1 [" a 3 I a 2 a s ] 
L k 3 I fci fc J L k 3 I k 2 k 3 J 

x I j J ao ai I a s \ [ a 3 I at ao ~| r a t I a 2 ai ] (D.44) 

— ~ / act t \ a 2 a 3 I a t J [ k 3 \ k t k J [ k t \ k 2 ki J ! 

S 

where S = ^ + zR+. The coefficients {. . . } can be represented explicitly by the following formula 



1 a a x I a 3 1 _ Sbjat + &o + «3 - Q)Sb{a3 + Qq - Qt)'S l b(Q3 + at - ap) 
Sb(a s + «2 + «3 - Q)Sb(aa + a 2 - a s )Sb(a3 + a s - a 2 ) 
Sb{a t + a 2 - ai) 



(DAS) 



S b (a s + a - ai 



•*4(U U U 2 , U 3 , U A ; Vi,V 2 , V3, V A ; 0) 



where 



U± = a s + ao — ai Vi = Q — a t + a s + 0:3 — oti 

U 2 = a s + Q — a — a>i V 2 = a t + a s + a 3 — ai 

U3 = a s + a 2 + 0:3 - Q V 3 = 2a s 

Ui = a s + Q!3 - a 2 V4 = Q. 

For completeness let us note that the b-Racah-Wigner coefficients { • • ■ } are related to the corre- 
sponding objects from 1IPT21 via 

v a s v a 3 

( olq a± a s 1 _ aittp a 2 a B ( a ai \ a B \ /p. in 

\ a 2 a 3 I a t / here — at „a 3 \ a 2 a 3 | a t / pT2 ■ (Ll.tO) 

L/ a 2 a\ ^atao 
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The b-Clebsch-Gordan and b-Racah-Wigner coefficients are meromorphic functions of all of their 
arguments. The complete set of poles may be described as follows HPT2l fBTp| 

{ a° al | a' } haf * Poles at 



— a 



32 s = — nb — mb 1 , a 3t0 = — n ^ — m ^ 



Q — ctgiQ = —nb — mb 1 , a[ 2 \ = —nb — mb 



l = 0,1,2,3, (D.47) 



where n, m G Z-°. 
[fc 3 3 | kt fe'] has poles at 



Q — = — n& — mb 1 , i = 0, 1, 2, 3, 

iife^cti+nft + mb" 1 , i=l,2, (D.48) 
± iks = Q — «3 + nb + mb -1 , 



where again n, m E Z- 



D.5 Proof of Proposition U 

Our proof of Proposition [4] will be based on the following nontrivial identity satisfied by the 
b-Clebsch-Gordan kernel. 



Lemma 10. 



fci ib JJ L 3 I fe 2 fci+ibj 



1 (D 49) 



T=-l 

where T^ 3 /(q?3) = /(as + t6), &3 = &2 + ki and we have furthermore used the notation 

&J=27ri Res [«s|E2], 

fcl — — IQi 



F r [£ 3]=27ri Res 



We have the explicit formulae 









ai 


Q!3 


-6 


«2 


ai 


Oi3 


-6 



_ e -*k 2ai *b\w2 — (D50 ) 



Sb{a 3 + a 2 -ai) Sb(ai + a 3 - ik 2 ) 



r Sb{2a 3 - 2b - Q) . 

^-La 3 -b = — — Tx — sin7r6(Q!2 + ai - a 3 ) , (D.51) 

^+ [ Z -I ] = S o { fo 3 ~JS sin T6(a 3 + aa + a x - Q) (D.52) 
x siri7r&(a3 + «i — 0:2) sin7rfo(o!3 + a 2 — ai) , 



7 We take the opportunity to correct some typos in these references. 
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where A a = a(Q — a). 

Given that Lemma [TOl holds, it becomes easy to complete the proof of Proposition [4] as follows: 
Notice that the left hand side of (ID. 491 ) can be written as 

Uai I 011 ~ b l\ la* I ° 2 ai 1 - r^bkt ["1 ~ ~ b )b r I a 2 ai ] 
LL 1 I fci »bJJ L 3 I fe2 fei+'iM - e 5 b (2ai +b) L I fc 2 fci+ibj 



e Trbk 3 



(K Sz ® E,J* • ["3 



ft 2 Ctl 1 
&2 fel J > 



where 0* denotes the transpose of an operator on 7^ ® 7^ defined by 

/" dkdk /(fcz,fcx)(Off)(fc z ,fc x ) = / rffc 2 dfc x (0*/)(fc a ,fc x ) ff (*; z ,fc x ). (D.53) 
Equation dP.49t may therefore be written in the form 

(K Si ® E Si )< • [a, I - £] = ^L SiSl (fc 3 )T; 3 [«3 1 - £], (D.54) 



where 

(kz) = e-^S b {2 ai + b)F T [T 3 %} [[«3 | 76 1 ■ 

By using the explicit expressions ( ID.5U .( |D.52I ) one may easily verify that the coefficients A^ 6 / 81 
are symmetric under the exchange of s 2 and s x , as claimed. This completes the proof of all the 
relevant statements of Proposition|4]for the case of the operator 02- 

In order to cover the remaining cases let us observe that 

2 = K S2 ® E 4l = B SiSi Oi B- 1 ^ , 

3 = R" 1 ^ (K- 1 ® E, J R S2Sl = B" 1 ^ Ox B S2Si , 

where B = PR is the braiding operator. 

Therefore, invoking Proposition [3] we conclude that Oi and O3 also satisfy Proposition [4] with 
coefficients A^if 2Sl (fc), r = 1,2,3 being related by 

A^(fc) = ^Sl^^fc) , = ^|-^^-(fc) , (D.56) 



respectively. The proof of Proposition|4]is complete. □ 

Proof of Lemma \TU\ Our starting point is the defining relation for the b-Racah-Wigner symbols, 
equation dP.44| i. Our claim will follow from (|D.44t as an identity satisfied by the residues of the 
meromorphic continuation of ( 1D.441 >. We need to analyze the relevant limits step by step. 

U\ = a s + ao — ol\ — ► —6: 
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Note that in the limit U± = a s + «o — &i —> —b the contour of integration in the definition of $ 4 , 
eqn. JA.14T >. gets pinched between the poles of the integrand s = Q — V4 = and s = — U\ — b as 
well as between s = Q — Vi + b = b and s = —U±. This implies that $ 4 has a pole when U\ = —b. 
In order extract the part which gets singular in the limit under consideration one may deform the 
contour of integration in ( IA. 14b to the sum of two circles around s = and s = b plus a contour 
which passes to the right of the pole at s = b and which approaches the imaginary axis at infinity. 
The residue is given as 

J_ . 2 S b (U 2 )S b (U 3 )S b (U 4 ) f smTrbUisimrbUzsmTrbUA 

27r Sm7r Sb (Vi)S b (V 2 )Sb(V 3 ) V sin 7r6Vi siii 7r6V a sin 7r6y 8 / ' ' 

Considering the behavior of {• • • } at U\ = —b, one finds that the pole of $ 4 at Ui = —b is canceled 
by the zero from the prefactor (S b (a s + ao — ai))^ 1 . Taken together one obtains the following 
special value for {• • • } at a s + ao — a.\ = — b : 

t a a 1 \ aa-ai-b \ _ S b (Q - 2a ~ b) S b {a t + OL 2 ~ OL X ) 



S b (2(ai - ao — b)) S b (a 3 + a 2 - a± + a + b) 

S b (a 3 + a t — op) S b (a 3 + a - a t ) 

S b (a 3 +a t — a — b) S b (Q + a 3 - a t - a - b) 

sin nbU 2 sin nbU 3 sin irbUi 



(D.58) 

+ a 3 - a t - a - o) 

( ein TrhTTn «in nrhT Qin ttJiTTa \ 

S b (a 3 + a + a t - Q) 1+ . . . . . 

\ sin7roVi smirbv 2 sm%bv 3 J 

The parameters U 2l U 3 ,Ui and Vi , V 2 , V 3 are now given by 

U 2 = Q - 2a Q - b Vi — Q - a t + a 3 - a - b 

U 3 = a\ — a — b + a 3 + a 2 - Q V 2 = a t + a 3 — a - b 
U4 = a\ — ao — b + a 3 — a 2 V 3 = 2(a 1 — a Q — b). 

k — > ia : 

In the same way as in the previous paragraph one may show that the b-Clebsch-Gordan coefficients 
[ k". I fci 1 fco ] and [ kl I kt to ] develop poles, with residues given by ( ID.50b . 
Continuation to Re(«o) = —b with Re(a3) = % — 6, < 5 < b < b^ 1 . 

The left hand side of JD.44b is analytic in the range under consideration. In order to describe the 
analytic continuation of the right hand side let us note that in the continuation from Re(ao) = Q/2 
to Re(ao) = — b exactly three poles at = a\ k \ k = —1,0, 1 cross the contour of integration, 
namely 



k = 


-1 : 


(fe) 
a\ 


= a 3 


+ a 


k 


= : 


(As) 
"t 


= "3, 




k 


= 1 : 


(As) 


= a 3 


- a a 



The analytic continuation of the right hand side of ( ID.441 I may therefore be represented by replacing 
the integration contour § in (ID.44b by C = § U Uk=-i ^k> with Ct being a small circle around the 
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i * ( fc ) 
poles at at = a t . 

Limit a,\ — > -6 with Re(a 3 ) = Q - 5, < 5 <b < b^ 1 . 

We observe that the integral over § vanishes due to the factor Sb{Q — 2a% — b). This is not the case 
for contributions from the poles at — a[ k \ k = — 1, 0, 1. Our claim now follows by straightforward 
computations. □ 



E. Construction of the Q-operator Q(u) 
E.l Preliminaries 

Let us now enter into the construction of the Q-operators. We begin by collecting some useful 
preliminaries. We will work in the Schrodinger representation where the operators x r , r = 1, . . . , N 
are diagonal. We will need operators U, Q., and J s defined in the Schrodinger representation by the 
following integral kernels 

N N 

C/(x,x') = n S(x r+1 -x' r ), fi(x,x') = n 5(x r +x' r ), (E.l) 

r— 1 r— 1 

N 

J s (x,x') = ( Wb ( l §-2 S )) N J] D s _ iQ (x r -x' r ). (E.2) 

r=l 

U is the cyclic shift operator defined in d82l . Q and J s are products of local operators, 

N N 

Q = l[n r , J s = l[}r. (E.3) 

r— 1 r— 1 

Here j r is the operator which intertwines at the site r the representations V s and V- s of 
Wq(s[(2,M)) (see 1PT21 '). and n r is the operator which realizes at the site r the parity opera- 
tion: ft r f(x 1 ,x z , . . . ,x r , . . . , xn) — f(xi,x z , . . . , —x r , ■ ■ ■ , xn) (whence n r Pr&r = 0(p r ) and 
n r x r Q r = 0(x r ), where the automorphism 6 is defined by d36ll). 

We will denote the standard scalar product on L 2 (M.) by (f\g) — L dx/(x) g(x). For a given 
operator 0, its transposed 0* and hermitian-conjugated 0* are defined, respectively, by (the bar 
denotes complex conjugation) 

(O'/ls) - (f\0g) , (0*713) - (7\0g) . (E.4) 
This definition extends to a matrix with operator-valued coefficients as follows 

(L%=(L V )\ {L% = (L v y t (E.5) 

i.e., component-wise. If is represented by the integral kernel 0(x, x'), then the kernels of its 
transposed and hermitian-conjugated are given by 

O t (x,x') = 0(x\x), 0*(x,x') = C>(x',x). (E.6) 
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In particular, we have 

U*=U* = U" 1 , ft* = ft* = ft, (E.7) 

4 = J s , j: = J-. = j; 1 . (e.8) 

In the Schrodinger representation we have x* = x, p* = — p, x* = x, p* = p and hence (as seen from 
©): 

4=e_ S) f* = f- s , ki = k:J, (E.9) 

Properties of the transfer-matrices of models in question with respect to the transposition and 
hermitian-conjugation are described by the following statement. 

Lemma 11. For the operations defined by AE.4\l and (IE.5\) we have 

( T ? xz ( u ))* = (-1) N T™(- U ) , (Tf (u)) 1 = T^i-u) (E.10) 

(Tf"^))*=Tf^), (Tf («))* =Tf («), (E.11) 

where T_ s (u) = J s T s (it) J~. 

Proof. Taking into account ( |E.9t , we observe that the L-matrices ( TTOb and ( f33l > satisfy 

(.-^s ( U )J =-cr 3 e L_ s (-u) e a 3 , 

(iS»)* = ff3 i a - .(-«)^ s 
Substitution of these relation into 

T*(u) - tr ((Ll)' • (L\r ■■■■■ (I&) T ) = tr (L*, L| • L\) (E.13) 

yields (IE. 10b . Relations jE. lib are derived analogously by noticing that we have (L(h)J = 
(73 L(m) (73 for both models in question. □ 

A consequence of this Lemma is that it suffices to prove Theorem |2] only for Q^(u). Indeed, using 
( 1E.61 > and (1A.241 >. it is easy to conclude that 

Ql(u)= (Z>_,(u)) N (Q>(u))*. (E.14) 

Therefore relations (|55]-i)-(|55]-iii) for Q_(tt) then follow immediately if we take ( IE. lib into ac- 
count. To check the Baxter equation d55Tiv) for Q_(m), we take hermitian-conjugation of (fSBTiv) 
for Q|(w), using (l55l-ii) and the property (IE. 1 Ob . After replacement of u by u this yields for 
Ql(it) = (Q+(it))* the following equation 



(E12) 



T b (u) • Qi(u) = (a(u))" Ql(u + »6) + (d(u))™ Q^(u - i&) . (E.15) 
Using relations dA.2 1 b and ( IA.23I ). we observe that 

D- s (u + ib) ——_ D- S (u- ib) 



a O) = — n — F^T~ d(«) = — n / \ "(")■ (E.16) 

D- S (u) D- S (u) 

Whence we conclude that (u) defined by (IE. 1 41 > satisfies (TSBTiv) (with the same coefficients a(u), 
d(u) as Q+(u) does). 
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E.2 Construction of Q-operators 

In order to construct the Q-operators explicitly let us consider the following general ansatz for the 
Q-operator: 

Q(u) = Y(«)-Z. (E.17) 

We will prove Theorem [2] for Ql(tt) in three steps: first constructing a suitable solution for Y (u) 
by requiring the Baxter equation to hold, then determining the form of Z , and finally checking that 
the obtained Q-operator satisfies (l55l . The first step in this proof is based on the idea to find such 
a gauge transformation of the L-matrix that it becomes effectively upper-triangular. This approach 
was originally applied by Pasquier and Gaudin MPG1 to the Toda chain. Our computation has many 
similarities with the modification of this approach developed in [De DKM| for the non-compact 
XXX magnet. 

Proposition 5. Let T b (u), b = xxz, sg be the transfer— matrices corresponding to the L— matrices 
( 1701 ) and A33t . Let Y b (u) be defined in the Schrodinger representation by the kernel 

N 

r„(x,x') = Y[ Di (u _ a) {x r -e h x' r+1 )D_ i (tt+0 .)(a; r -<), (E.18) 

r=l 

where £ xxz = 1, s SG — — 1. Then Q b (u) of the form ( IE. 7 71 ) satisfies the Baxter equation ( I55l -iv ) 
with coefficients a(u), d(u) as specified in eq. (162 P of Theorem]2\ 

Proof. Let us introduce the gauge-transformed Lax operators (the transformation depends on the 
site number r): 

L b r (u) = G b r+1 ■ L b r (u) ■ (Gly 1 , Gt=Q?), (E.19) 

p xxz = e nb(2 X ' r -u) ( p SG = e 2,b x ' r ( g 2Q) 

The relevant matrix elements of the new Lax matrices are given by 

(£r(«)) 21 =4*b(/>rPr+i)* (E.21) 
x ^cosh7r6(x r — e v x' r+l + ~(<r — u)) cosh7r6(x r — x' r + |(<t + u)) k r 

- cosli7r&(x r — e\, x' r+1 — \{o — u)) cosh"7r&(x r — x' r — \{a + u)) k^ 1 ^ , 
(2£(u)) =2^e^- £bX " ) ( e **nr&(t»-«r) cosh7r6(x,. -4 + i(cr + M))k,. 

- e$ e > *'(»-«0 cosh 7rfe(x r - x' r ~ \ (a + u)) k" 1 ) , (E.22) 
{pM)) 22 = 2^e^<+^ (e^ 6 ^) coshTrbixr-ebx'^ + ^u-a))^ 1 

- e ~^ b{u+a) cosh7r&(x r - e b < +1 + \(a - u)) k r ) , (E.23) 
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where x xxz = 1, k sg = —ie . 

-l 



In the Schrodinger representation, operators k r , k r act as shifts of x r by ±|6. Using the functional 
relation ( IA.2 11 >. it is straightforward to apply ( 1E.2U to Y„ (x, x') and verify that the condition 

(Z^(u)) 21 F tl b (x,x')=0 (E.24) 

is satisfied for all x' G R N . This implies that Z b (u) becomes upper triangular when acting on 
Y b (x, x') so that we can calculate the action of M b (u) on Y b (x, x') as 

M b ( W )^(x,x')- fn^i(4(«)) u * W(x,x'). (E.25) 

Hence, taking into account the periodicity condition, G^ +1 = G b , we have 

(N N \ 

n( z ^))ii+n( z ^)) 22 Kwo- (e.26) 
r=l r=l / 

Applying (IE221 . (IK231 to Y u b (x, x') and using dE.241 >. we derive 

(Z b (u)) n F u b (x,x') = 2e b e ^<-< + ^ S inhnb(u - a) y„ b _ i6 (x, x') , 
(Z b (u)) 22 y u b (x,x') =2 4 e^(^-<' sinh7r&( W + a)y u V i6 (x,x'). 



(E.27) 



'22 

Here we have used that Y^(x, x') = Yl r =i ^( u > x r)> where each factor satisfies the relation 

Y b (u,x r ±i^) f cosh7r&(a; r — x' r — + u)) \±i 



r±!^J / COSH 7tD(Q 

Y b (u ± ib, x r ) vcosh7r&(x r — £ b x' r +i + 2"( CT — u )) 
Combining ( IE. 26b with ( IE.27I ). we obtain 



(E.28) 



T b ( U )Yi(x,x') = {2x b S inh 7 r&( U + a)) N r„ b +4b (x,x') 



+ (2e b xr b sinh7r&(w - cr)) N Y b _ tb (x, x') 



(E.29) 



which implies that the Baxter equation d55l-iv) holds with the coefficients a(u), d(u) as specified in 
eq. d62j>. □ 

The possible form of Z can be found from the requirement that (l55l-iii) holds. 

Proposition 6. Let Y (u) be chosen as in Proposition^ Then the commutativity condition 

Q b (u)T b (u) =T b (u)Q b (u) (E.30) 

holds for Q b (u) of the form iE. 1 71 ) provided that the corresponding operator Z b satisfies the following 
relation 

Z b T b (w) = T b _ s (w)Z b . (E.31) 
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Proof. In order to treat both models in a uniform way, let us introduce the operator 

i 1 , \> — xxz 

Q b = { , (E.32) 

where the parity operator was defined in ( IE. 11 1. Then, using the explicit expressions dE.18b and 
taking ( IA.22b into account, it is easy to verify that (the subscript x or x' of an operator specifies the 
argument on which it acts) 



fi x *2(x, x') = q x , y„ b (x, x') = u- 1 y^(x', x) . (E.33) 



Now we derive 



T s b ;x ( u )rj(x,x') ^3 ( a(M )) N r^ b (x,x') + (d(«)) N C(x,x') 



jE33l 



(- £b ) N ^ u x 1 ((d(-u)) N ^ u+i6 (x', x) + (a(- M )) N yl u _ i6 (x', x)) 
£b ) N Yl (x, x') (T- x , (-«)) 4 ^ (x, x') Tl s;x , («) . 



Thus, we verified that T b s (u) Y(it) — Y(w) T b _ s (u), which is equivalent to (IE. 30b if relation (IE. 31b 
is satisfied. □ 

Proposition|6]implies that we can choose 

Z b = J s , (E.34) 
where J s was defined in dE.2l ). Then we compute the integral kernel of Q^_(u): 

r N 

Q+. u (x,x') = (wb{i® - 2s)) / dz\ . . .dz N D^{z r - £\,x r -i) (E.35) 

1 — i 

x Z3_ (z r — x r ) D s {z r — x' r ) 

N 

= (D- S (u)) Ds^(x r - x' r ) Dv±u(x r -i - e b x' r ) D- S (x r - e b x r -{) ■ (E.36) 



Equivalence of ( IE. 35b and ( IE. 36b is due to the identity dA.341 >. 

Remark 9. The Baxter equation d55l-iv) along with the self-commutativity (|55l-ii) relation, which 
will be proven below, imply that (IE. 30b extends to commutativity of Q+(it) with T b (v) for those 
values of v, where Q+(w) is invertible. 
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E.3 Proof of commutativity relations 

To complete the proof of Theorem[2] we have to establish relations (l55l-i) and OBl-ii). 
Lemma 12. Let Y (it) be chosen as in Proposition^ Then the following identities hold 

(Y b («))*. Y b (v) = (D s (u) !>_,(«)) N (Y b («)) * • Y b («) , (E.37) 
Y b ( u )- (Y b («))* = (£>,(«)£_.(«)) *V(i;)- (Y b (S))\ (E.38) 

Proof. These identities are just particular cases of the integral identity dA.351 . Indeed, let us denote 
a u = h(u— a), j3 u = —\(u + a). We will also use the notation a* = —\Q — a. Let us consider 
the operator V b (it, v) — (Y b (it)) * • Y b (v). Its kernel is given by 

N 

T/ b U!t ,(x,x') = / dzi...dz N Y[D a *(z r — ebX r+1 )Dp*(z r -x r ) (E.39) 



R N 



r=l 



xD av (z r - e b x' r+1 ) Dp v (z r - x' r ) , 



Now we can apply identity dA.35b choosing a = a*, (3 = a v , 7 = /?*, w = f3 v , and u = e\, av+i, 
v = e\, x' r+1 , w — x r , z — x' r . This yields 



N 



^(x,x')=(A(a>„,^A)) (e b (x r+1 -x' r+1 ))D^(x r -x' r ) 

r=l 

N 

x dzi... dz N Y[ D au (z r - e b x' r+l ) Dp u (z r - x' r ) 



r=l 

IK."" 

X _D Q * (z r — E\j Xr+i) Dp* (z r — Xr) 

= {D s (u)D^ s (v)f V v \ u (x, x') . (E.40) 

Here we used (1A.221 I. the definition ( IA.28b of the function A(a±, a.2, ■ ■ .), and took into account the 
periodic boundary conditions. 

Identity ( |E.38t can be proven absolutely analogously. □ 

Proposition 7. The operators Qj_(it) and Qt_{u) with the kernels given in Theorem\2\by eqs. ( 1601 ) 
and ( 1671 ), respectively, satisfy the following commutativity and exchange relations 

Q b + (u)Q b + (i;) = Q b + («)Q b + (u), Q b _(u) Q b _(i>) = Ql(v) Q b _(u) , (E.41) 
Q b + (u) Q b _ (v) = Q b _ (u) Q b + (v) = Q b + (v) Q b _ («) = Q b _(v)Q+ (u) . (E.42) 



for all u, V G C. 
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Proof. Observe that, using ( IE.6b and ( 1A.241 >. the equality of (IE. 35b and ( IE. 36b can be written in the 
following operator form: 

Q\{u) = y\u) ■ J s = {D_ s (u)f (X^ 1 • & • IT 1 ■ (Y\u))* , (E.43) 

where we used the notation ( IE.32I ) and introduced 

N 

X^ = JjD^Xr-efcXr-i). (E.44) 

r=l 

Using Lemma[T2l we can write down the product of two such Q-operators as follows 

Q b + (u) Q» + (v) = {r>-s(u)f (X^ 1 • Q b ■ IT 1 • (Y b (u)Y ■ Y b (v) ■ J s 

(D_ s (v)) N (Ky 1 ■ fi b • IT 1 • (Y\v)Y ■ Y b (u) ■ j. = Q b + (v) Q>) . 



(E.45) 



This proves the first relation in dE.41b and hence dB31-ii) for Q\(u). As was explained in Subsec- 
tion lE.ll relation d33]-ii) for (u) (i.e., the second relation in (IE. 41b ) follows then as a consequence 
of the relation ( IE. 14b between Q + (u) and Ql(u). By the same token, relation (IB31-i) is equivalent 
to ( IE. 42b . To prove the latter relation, we substitute ( IE.43b into (IE. 14b and use ( IE.7b - (lE.8b . This 
yields the operator Q b _(u) in the following form: 

Q b _(u) = {D_ s (u)) N J; 1 • {Y b (u))* = Y b (u) ■ U • fi b • Xt . (E.46) 

As seen from (IE.43b and (IE.46b . the two expressions on the l.h.s. of (IE.42I ) are just two ways to write 
down (D_ s (i>)) Y b (u) ■ (Y b (w))*. Analogously, the two expressions on the r.h.s. of ( IE. 42b are 

two ways to write down (D_ s (ii)) N Y b (v) ■ (Y b (uj)* , The middle equality in ( IE. 42b is due to the 
identity dE.38b in LemmaUll □ 
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